For Reference 


NOT TO BE TAKEN FROM THIS ROOM 


Gx ipnis 
UNIMASITATIS 
HARERTAEDSIS 


Digitized by the Internet Archive 
in 2024 with funding from 
University of Alberta Library 


httos://archive.org/details/Olaosebikan1976 


we bi 
13 ha 


bara 


i. seyeliak 


AS ihe ney 


ae 
Ea ih ae 


THE UNIVERSITY OF ALBERTA 
RELEASE FORM 


NAME OF AUTHOR : OLALEKAN JACOB OLAOSEBIKAN 
PET oR OF Sloe oho s ON THE VIBRATION OF ELASTIC 
SPHERICAL SHELLS | 
DEGREE FOR WHICH THESIS WAS PRESENTED: MASTER OF SCIENCE 


YEAR THIS DEGREE GRANTED: 1976 


Permission is hereby granted to THE UNIVERSITY 
OF ALBERTA LIBRARY to reproduce single copies for 
private, scholarly or scientific research purposes 
only. 

The author reserves other publication rights, 
and neither the thesis nor extensive extracts from 
it may be printed or otherwise reproduced without 


the author's written permission. 


eas 
an ; 


egy “uakrosne aqanad — 


aS ry sisi 2 ; 4 Fi qui ad és fab bade ; 
ity; ca iva Wer. 
ais ce i Pearl ae ie 2 63 —enra).’ ATA : TO y 


di) Goagpti< oftirueley Ge eiatogan’ meters 


THE UNIVERSITY OF ALBERTA 


ON THE VIBRATION OF ELASTIC SPHERICAL SHELLS 


by - 


OLALEKAN JACOB OLAOSEBIKAN 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF MASTER OF SCIENCE 


DEPARTMENT OF MECHANICAL ENGINEERING 
EDMONTON, ALBERTA 


SPRING, 1976 


a | 
a 


nized: ss ae nisaiinaatts a 


7. 


7 f.' 
igagirs Ife Auf nos, raf des oe | ~ 
, 0 / ¢ 


ie 
pi cany A 
j s 
, nh = 
PAT Cid, PTS Hitocdyd to Yoauadd or qe 
amon Ne ner aula Sir PRMAee SATAT OE 


Sow Tos (a0 ur eT) 


; : QvGtanie o4-clAngen qo TeaNY festa 
| ROMA KOT eT 
4 APR coins 


THE UNIVERSITY OF ALBERTA 


FACULTY OF GRADUATE STUDIES AND RESEARCH 


The undersigned certify that they have read, and 
recommend to the Faculty of Graduate Studies and Research, 
for acceptance, a thesis entitled "On the Vibration of 
Elastic Spherical Shells" submitted by Olalekan Jacob 
Olaosebikan in partial fulfilment of the requirements for 


the degree of Master of Science. 


, in.) pil é 
» avert yeah decty Aries bonbbaidieay weit | a c 
ieesek hae bate «teu % stioeet ale 8 ieee 7 
te im Lawn? pdt AN” Laer da gevdike 3 deumnagetom 208 a. ; 
den win hele ya coll gpa ate eh inte Lil Migalig® orseate * a 
Pra | ERPS ETN pig T° heed Pee 472 Swlowaey oe metideeneio nak 
Se i ATT e pnt wh Bi 
; > ss 


- — 


ABSTRACT 


The transient response, to time dependent pressure 
and velocity inputs, of elastic spherical shells made of 
an isotropic and homogeneous material is investigated. 
Standard numerical methods are used to solve the governing 
equations for the incompressible thin and thick shell and 
for the compressible thin shell problems of non-linear 
elasticity while a numerical procedure is developed for 
the method of characteristics to solve the compressible 
thick shell problem of linear elasticity. The results 
obtained for the free and forced vibrations of the thin 
and thick shells considered agree well with existing 


analytical and numerical solutions. 
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LIST OF SYMBOLS (continued) 


2. Nomenclature for the Thick Shell Problem: 


Incompressible Case of non-Linear Elasticity 


(A,a) 


(Bb) 


(R,r) 


=I 
rl 
Fle 


inner radius in the (undeformed, deformed) 
state 


outer radius in the (undeformed, deformed) 
state 


radial coordinate of particle in the 
(undeformed, deformed) state 


radial displacement at reference radius, r 
time 

io yi ; : 

z = non dimensionalized time 


density 


shear modulus for infinitesimal strain from 
the undeformed state 


strain energy function 


non-dimensionalized strain energy function 
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LIST OF SYMBOLS (continued) 


Nomenclature for the Thick Shell Problem: 


Compressible Case in Linear Elasticity 


c = wave speed 
) = density 
Aa = Lame's constants 
Vv = time 
1s = radial distance 
u = radial displacement 
ou he : : 
Mek Tare = particle velocity 
(A,B) = (inner, outer) radius 
an = radial stress 
4 = tangential stress 
o = 4 = = non-dimensionalized stress 
r= = = non-dimensionalized radial distance 
T = (B-A) = shell thickness 
GS z = non-dimensionalized particle velocity 
u = a = non-dimensionalized radial displacement 
T= oe = non-dimensionalized time 


E = 2(1+v)u 


Young's modulus 
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LIST OF SYMBOLS (continued) 


(l= eae 
2 GEVE) = dimensionless parameter 
o, - o0_ = discontinuity in stress, 6 
Vv, - V_ = discontinuity in velocity, v 


= value of stress, o, ahead of the wavefront 
= value of stress,.o, behind the wavefront 


= value of velocity, v, ahead of the wave- 
ELON 


= value of velocity, v, behind the wavefront 


XxX 


alee ohn mar 


i _s —— es =r Mealle Pa 


nn yee 


sma NE yo ad aa 
ae) has te hann “ Seon 


sasdnatl wade Wat igual © 4 ine wo qin 3 i : 


CHAPTER 1 
INTRODUCTION 


In the analysis of the vibration of elastic 
spherical shells of isotropic and homogeneous materials, 
‘different cases of the dynamical problem arise according 
to (i) whether the linear or the non-linear theory of 
elasticity is considered; (ii) whether the material of the 
shell is compressible or incompressible and (iii) whether 
the shell considered is thin or thick. Free and forced 
vibrations are usually treated as two aspects of the same 
problem resulting from different initial and boundary 
conditions. 

The methods which have so far been used to solve 


these problems are outlined for the different cases as 


follows. 
1. Free oscillations of linear elasticity theory 


The linear elasticity problem of the free rad- 
ial oscillation of the incompressible and compressible 
thick shells, with the thin shells considered as lim- 
iting cases, have been solved by means of spherical 
harmonics. This method of solution is discussed fully 


by LOVE [1]. 
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Forced vibration of a compressible shell in linear 


elasticity theory 


The forced vibration problem of linear elast- 


icity for a compressible thick shell, again with the 


thin shell considered as a limiting case, has been 


solved by a number of analytical and numerical methods. 


The analytical methods of solution include 


(i) 


CEv) 


Caa35) 


(i) 


(iy) 


Gra) 


(iv) 


a direct infinite series solution using the 
Midlin-Goodman technique (BAKER and ALLEN [2]), 
a Laplace transform method (TRANTER [3]), and 


a finite Hankel transform method (CINELLI [4]). 
Numerical methods of solution ines 

the discontinuous step method (MEHTA and 
DAVIDS [5]), 

a finite difference method (SMITH [6]), 

a combined characteristic-difference method 
(CHOU and GREIF [7]), and 

the method of characteristics (CHOU and 


KOENIG [8]; LEONARD and BUDIANSKY [9]; 


CHOU and GREIF [7]). 
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Forced vibration of the incompressible thin and thick 
These non-linear elasticity problems of the 

thin and thick incompressible shells have been solved 
for admissible strain energy functions (ERINGEN [10], 
ZHONG-HENG and SOLECKI [11]). A well-known admissible 
strain energy function for an incompressible material 
is the Mooney-Rivlin type [given by equation (2.6) of 
Section (2.2)] from which the neo-Hookean form is 


obtainable as a special case. 


Forced vibration of the ythin compressibiesshel lyot 
non-linear elasticity 


This problem has also been solved for ened 
Sible strain energy functions. The strain energy 
function considered here is the Blatz-Ko type [given 
by equation (2.8) of Section (2.2)] proposed for a 
Class of rubbery materials. It gives the neo-Hookean 
form as a limiting case when the value of Poisson's, 
ratio, v, approaches 0.5. A number of other strain 
energy functions, which are admissible for a compres- 
sible material, have been discussed by HADDOW and 
FAULKNER [12]. 

In this investigation, the following free and 


forced vibration problems are considered. 
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(i) The incompressible thin shell problem of non- 
linear elasticity theory for the Mooney-Rivlin 
and a logarithmic strain energy functions. 

(ii) The compressible thin shell problem of non- 
linear elasticity for the Blatz-Ko strain energy 
function. 

(iii) The incompressible thick shell problem of non- 
linear elasticity theory for the neo-Hookean 
solid. 

(iv) The compressible thick shell problem of linear 


elasticity solved by the method of characteris- 


tics. 


For forced vibrations, a step input applica- 
tion of pressure is considered but only a little mod- 
ification of the governing equations is required for 
other types of pressure and velocity inputs. 

For all the problems considered the governing 
differential equations of motion, expressed in non- 
dimensional form for generalization, are integrated 
by suitable numerical methods which are discussed in 
the appropriate sections of this thesis. 

In order to verify the results of these anal- 
yses, various comparisons are made between the results. 


For example, the results for the period of free radial 
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oscillation of the neo-Hookean thin shell, obtained 
separately from the incompressible thin shell theory, 
the compressible thin shell theory and as a limiting 
case of the Hookean thick shell theory, are compared. 
The computed results are also compared with the sol- 
utions obtained by other investigators for a number 
of different boundary and initial conditions of the 
forced oscillation of the thin and thick shells. 

In all these comparisons, the results, shown 
in graphical and tabular forms, are found to be con- 
sistent. 

Finally, for completeness of the investigation, 
the variation with time of the tangential and shear 
stresses at the inner surface of the compressible 


thick-walled shell are plotted as shown in Figures 20 


and 22. 
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CHAPTER 2 


THE HYPERELASTIC MATERIAL 


Strain Energy Functions and Constitutive Relations 


The strain energy function, W per unit unde- 
formed volume, of an isotropic hyperelastic solid 


is a function 


W = W(I,, ig T,) (ed?) 


By 


of the basic invariants 


| 
II 
cr 
i) 
w 
H 
II 


$[ (tr B) 2 - tr B?] and 


I, = det B (242) 


of the left Cauchy-Green strain tensor B = FF 
where F is the deformation gradient tensor. 


The constitutive equation is 
2 
¢=pi+ (6 + I,~)B- pB os) 


where o is the Cauchy stress tensor and 


6 OW -'5 OW -; 9W 
p=21,7 =, ¢=21 a=, V=2I raped 
3 sr | dr, 3 dr, 

(2.4) 


The constitutive equation (2.3) reduces to 


the classical stress-strain relation 
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Of=—25llee + (4 tre). Lf (Pres) 
when terms area are neglected. 


When the material is incompressible, I. = ih 
and p becomes a Lagrangian multiplier, which is no 
longer determined by the deformation but must be 


obtained from equilibrium conditions. 


2.2 Particular Strain Energy Functions 


We consider the following strain energy fun- 


ctions 
W =$ Velen se Ube) Ge ay) anes thal 
(2.6) 
u ae 
Wee 5 [c(I,-3) + 3(1-c)2n 5 ] where C30 
(7) 
and 
ail Peay a SL 
Wes | (3) = bs i) 3 (258) 


where uw is the shear modulus for infinitesimal 


strain from the undeformed state and 


Ss = ay (2.9) 


where v is Poisson's ratio for infinitesimal 


deformation from the undeformed state. 
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The strain energy function given by equation 
(2.6) was proposed for an incompressible material 
independently by Mooney and Rivlin? Equation (2. /) 
is also for an incompressible material and was pro- 
posed jointly by Gent and Thomas , while the function 
given by equation (2.8) was proposed for a class of 
compressible rubbery materials jointly by Blatz and 
Ko’. This latter strain energy function has been 
examined in detail by HADDOW AND FAULKNER [12] for 
the static problem of finite expansion. 

A special case of each of these strain energy 


functions is the neo-Hookean form 


W = 5 (i so Sh)) (2rn 0) 


which is obtained from equations (2.6) and (2.7) when 
c= 1 and as a limiting case of equation (2.8) as v 
approaches 0.5. 

The restrictions on W, which must be satisfied 
to ensure a periodic motion of the spherical shell, 


are discussed in section (B.1) of Appendix B. 


* 
See ADKINS [15]. 


a See BLATZ and KO [16]. 
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CHAPTER 3 


THIN SHELL PROBLEMS 


Approximate Theory for Thin-walled Shells 


Using the notations for the thin shell problem, 


let 


(3.2) 


All 


and assume that € is small enough to justify the 
neglect of terms leas and this is done without com- 
ment in what follows. The mean radial and circumfer- 


ential stretches are then given by 


eee mand? mie oS (seo) 
T R 
respectively, and 
2 ees ee) 
iG ea 
The following relations, then, hold 
Be=eRe(L PhS), PARA =aRe (ieee) (3.4) 
and 
= ce} — 50) 
joy ee ie 1 Ge eee a=r(1l- —=) (sii) 
on ! 2X 


Mean radial and circumferential stresses are 


defined by the equation 
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a 
(3.6) 
This equation, along with equations (3.5), gives the 


relation 
S. = - 5 (a, +a) (337) 
2 ik Outs ; 


The equation of motion of the thin shell can 
now be derived from the theorem that relates surface 
tractions to the mean hydrostatic stress as follows. 


This theorem, the proof of which is given in 


Section D.1 of Appendix D, leads to the equation 


| T- yrds - | or * r dV = 3Vo (3.8) 
Cee V 


where V is the volume, S the surface, T the surface 


traction, r the position or radius vector and o is 


the average hydrostatic stress. 


It follows from this last equation, along with 


equations, (3.5)mand (8.6), that 


= il nN oa 
C= 5 [ = (q,-q,-ert) - (q,t+q,)] 
(3.9) 


and on substituting equation (3.7) we get 
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This is the equation of motion for the thin-walled 
shell of either compressible or incompressible mater- 
ial, the final form of which is obtained after sub- 
stituting for oO, and oe from equation (2.3) for a 
given strain energy function, W. 

For convenience, we drop the bars on the mean 


quantities so that equations (3.7) and (3.10) become 


respectively 
6) te £ (qua ct) (Sel) 
te 2 Hh re) c 
and 
eaeee a Ge g ipet) (3e1 2) 
9) ie 2E0 ny re) x 


where p = p, / (647) is the mean density of the shell, 
p, is the density in the undeformed state, and r is 
the second derivative of r with respect to dimension- 
ed time T. 


For a thin-walled spherical shell, the mean 


physical components of the tensor B are given by 


0 0 d . (S23) 
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The Incompressible Case 


Incompressibility condition gives 
6 =} (3551/4) 


Substituting equations (3.13) and (3.14) in equation 
(2.3) and putting the resulting expressions for oe 
and oir in equation (3.12) we obtain the equation of 


motion 


Spel Chg es = 2el(y = 76 +0 - <3) cE 


(35115) 
To non-dimensionalize this equation, we divide it 


tEheoughn by Wwato obtain 


q. - q ee 
Ze 0 prt 1_ 1) ~ t)¥ 
(3416) 
Ti 7 Wo : . ; 
Putting i = Q as the non-dimensional input 


pressure at the inner surface, 
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W ; A 
where r is the second derivative of r with respect 


to the non-dimensional time, T, then equation (3.16) 
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(3) 

On substituting values for % and ¥ for a given strain 
energy function W, subject to the incompressibility 
condition (3.14), equation (3.19) gives the final 
non-dimensional form of fie equation of motion of the 
thin incompressible shell. It is a second order non- 
linear ordinary differential equation which can be 
solved numerically by using for example, the fourth 
order Runge-Kutta process or a modified version of 
the method of discrete variables as described in 
Appendix A and shown in Programs 3a and 4a. 

We note further that after substituting for 
6, ¥ and Q, equation (3.19) can be integrated direct- 
ly to obtain the energy equation. Thus, for the 


Mooney-Rivlin material, the equation of motion is 
y" = QR _ 2 ein - 4) = 2(a-e) (A7- 4) 
€ Ne 3 


(3.20) 
£Or constant O. 


This equation gives 
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Integrating this last equation and using the initial 


conditions 
| 
MeO) = 0 and A(0) =1 iC Epe ae 


we get the energy equation 
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The conditions for the motion of the shell to be 
periodic follow at once as discussed in Section B.1l 
of Appendix B. 

In the special case of the neo-Hookean solid, 
for which c = 1 in the foregoing equations, equations 
(3.20) and (3.23) reduce respectively to 
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To find the maximum value, Ort! of the applied 


step input pressure Q, for which a periodic motion is 
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possible, we consider the critical case depicted by 

curve (b) of Fig. 2, for which hg is a common root 

of equations (3.24) and (3.25). Following the pro- 

cedure of Section B.3 of Appendix B, we solve these 
' t 


equations simultaneously, when 4 = 0 and i =0 


respectively, to obtain for e« = 0.05, Ce = 0.0556 


at ha = 1.7317. This critical condition is shown in 
Fig. 10b. Furthermore, for a value of Q < Ore 
say Q = 0.02, the finite period of oscillation is 
2.0864. The response curve for this case is shown 
in Fig. 14a. Also, as Q approaches zero, say Q = 10 
the period of oscillation is approximately 1.8141. 
This is the same as the result of linear elasticity 
theory as further discussed in Section 6.3. The 
response curve for this limiting case is shown in 
Fig. 14b, and the phase diagram in Fig. 10c. This 


latter curve is elliptic, which indicates a harmonic 


oscillation, as further discussed in Section 6.3. 


The Compressible Case 
Let Cie ae and 2 ee 0, then equations (3.11) 


and (3.12) become respectively 
ee 3.26 
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For the Blatz-Ko material, the strain energy function 


W, given in equation (2.8), is 
an bi _i = 
Wea C(I, 3) res 1)] 


where 


1=2v 
and v is Poisson's ratio for infinitesimal deform- 
ation from the undeformed state. From this expres- 


sion and equation (2.3), we have 


(3.28) 


and 


(37529) 


Q 
Il 
On| 
J 
~ 
for) 


Substituting equation (3.28) in (3.26), we get 
q = -21 2 - (078)78? | (3.30) 
r 


and substituting equations (3.28) and (3.29) in 


(3.27), we have 


i] 2 52 
r = i SES es (3.31) 


Equations (3.30) and (3.31) are a pair of non-linear 
ordinary differential equations, in \ and 5, govern- 


ing the motion of a thin-walled shell of 
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compressible Blatz-Ko material. This pair of equa- 
tions is integrated numerically by combining the 
Newton-Raphson root-finding method with the method of 
discrete variables as described in the algorithm of 
Section A.5 of Appendix A. 

The static curves of Fig. lc are obtained by 
solving equations (3.30) and (3.31) simultaneously 
LOT ‘ = 0. These curves are the same form as the 
static curve (i) of Fig. la for the neo-Hookean mat- 
erial. Hence the conditions for a periodic motion 
of the thin shell of compressible Blatz-Ko material, 
are the same as given for the neo-Hookean thin shell 
in Section B.1 of Appendix B. 

The neo-Hookean solid is a limiting case of 
the Blatz-Ko material as Poisson's ratio, v, approach- 
es 0.5. For v = 0.49 and € = 0.05, results computed 
from equations (3.30) and (3.31) agree well with 
those of Section 3.2. Thus, Q = 0.0556 and 


max 


dg = bo73a20. Also, «or Q)=)0.02, >the finite; period 


of oscillation is 2.0893 and as Q > 0, the period is 


approximately 1.8157. 
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Specific Examples 


Based on the foregoing analyses, the follow- 
ing specific examples have been evaluated. 
(1) The neo-Hookean solid, 
(ii) the Mooney-Rivlin material for which c = 0.9 
and 


(iii) the Blatz-Ko material for which v = 0.3. 


3.4.1 The neo-Hookean solid 


The neo-Hookean solid is a special case of the 
Mooney-Rivlin material and a limiting case of the 
Blatz-Ko material. Results of numerical calculations, 
based on the strain energy functions for each of 
these latter materials, agree well. Thus, as shown 
previously in Sections 3.2 and 3.3, One = 0.0556, 
hg = 1.732, the period of oscillation for Q = 0.02 


is approximately 2.088 and as Q + 0, the period is 


approximately 1.815. 


3.4.2 The Mooney-Rivlin material 

Consider the case c = 0.9 in equations (3.20) 
and (3.23) when « = 0.05. Computation, described 
in Program 2a, gives 2.4, = 0.0633. This is shown 


in the phase plane diagram of Fig. lla. For Q <Q... 
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say Q = 0.06, the phase plane trajectory, of Fig. llb, 
closes, so that the motion is periodic with a finite 
period of oscillation computed from Program 5a to be 


4.2455. 


3.4.3 The Blatz-Ko material 


Consider the case v = 0.3 in equations (3.30) 


and (3.31) when € = 0.05. Computation, described in 
Program 2c, gives Oras = 0.047545. This is shown in 
the phase plane diagram of Fig. 13b. For QO < Oh aed, 
say Q = 0.02, the phase plane trajectory of Fig. 13c, 
closes, and so the motion is periodic with a finite 
period of oscillation, 2.6925. Also shown in the 
phase plane diagram of Fig. 13b is the case Q = 


0.06 > Oy, for which the phase plane trajectory 


does not close and so the motion is non-periodic. 
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Free Oscillations 


For free oscillation, Ceara sas 0, and the 
initial conditions are, in general, 
q i] 
ECO) ep and X8(0) = \ ee (Siz) 
If the initial perturbation is a purely radial 
displacement or velocity, then we have the simpler 


boundary conditions 


A(0) = A, and d (0) = 0 (333) 
and 

A(0) = 1 and d* (0) = a (3.34) 
respectively. 


When the amplitude of oscillation approaches 
zero, the response approaches that predicted by 
linear elasticity theory. 

For the Mooney-Rivlin material, the equation 
of motion, for free oscillation, obtained from 
equation (3.20), is 


r = -20(d - 5) jy = yee a Ss). 
r 


(3335) 
This is a second order non-linear ordinary different- 
ial equations which can be integrated numerically as 
before, subject now to any of the initial conditions 


(3532) 3= (3534) 
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In particulan,, Lor ¢ = 1, and with the anitial 
conditions: 1. = 0 and X,= 1, say 1.005, equation 
(3.35) describes a free radial oscillation of the 
neo-Hookean thin shell. Computations yield the phase 
plane diagrams of Fig. 10d. This phase trajectory 
is an ellipse, which indicates that the motion is 
harmonic. The period of oscillation is 1.8140, which 
is the same as the result obtained from the linear 


theory*. 


* See Rayleigh's method of Section 4.2. 
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CHAPTER 4 
THICK SHELL PROBLEMS 


General Discussion 


The ratio ¢€, of wall thickness T, to the mean 
radius R, of a shell is given, in the notation for 


thick shells, by the equation 


B 
ae : : vo (4.1) 
A 
where A and B are respectively the inner and outer 
radii of the shell in the undeformed state. When ¢ 
is sufficiently large that terms O(e*) are no longer 
negligible, the spherical shell is considered thick- 


walled. This condition usually holds for shells for 


which 2 sso Wehner = el ye 


For such thick-walled shells, the following 


problems are considered. 


(i) Free radial oscillation of the thick incompres- 


sible shell of linear elasticity using 
Rayleigh's method. 

(ii) Forced oscillation, under a step function 
application of pressure, of the thick com- 
pressible shell of linear elasticity using 


the method of characteristics. 
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(iii) Finite oscillation of the thick incompressible 
shell of non-linear elasticity using a strain 


energy function and an acceleration potential. 


The following special cases are then examined 


in Section 4.5. 


(i) The thick shell results approaching those of 
the thin shell theory when : <707 i. 
(ii) The results of non-linear elasticity approach- 


ing those of linear elasticity as deformation 


approaches zero. 


Incompressible Shell in Linear Elasticity: 


Free Radial Oscillation by Rayleigh's Method 


Consider a shell of incompressible material 
with inner and outer radii, A and B respectively, 
as shown in Fig. 3. For radial oscillation, this 
body is a single degree of freedom system and since 
there is no dissipation of energy, the single or 
fundamental frequency of vibration can be determined, 


using Rayleigh's method, from the energy equation 
uk = V (4.2) 


where LN er is the maximum kinetic energy and Ves 


is the maximum potential energy of the system. 
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To find Th , let a be the amplitude of 


ax 
oscillation at the inner radius A, and u_ the ampli- 
tude at the reference radius R, then, since the motion 
is harmonic, the maximum velocity at R is (uw), where 
Ww is the fundamental frequency of free oscillation 

in radians per second. If dm is the mass of the 


elemental shell of radius R, thickness dR, and 


volume dV, then 


B 
- 1 2 
Tae = | OY (uw) am, Cae) 
A 
and 
si = ei = Cnereain (4.4) 


where 0 is the density. It follows then from. 


equations (4.3) and (4.4) that 


ie oe 


B 
Te = [ 450u u R'dR. (425) 


For an incompressible material, a and u are 


related by the equation 


4 


4 n(n? - a?) =F m[(R + a)? - (A+ a) 3] 
(4.6) 
Neglecting terms Otacnn this equation gives 
MEAs ; (4.7) 
R 


Substituting for u from equation (4.7) in (4.5), 


we get 
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2 Ae? dR 

T = 2Tpw Aa | — 

max A R2 
Thus 

- 2 Ave. 1 il 
ie ee = 2T0w Aa (- B + x) : (ArG) 
TOs na Mieco we consider the strain energy V, 
per unit volume given by 
pag ; 
Nias on 54 (4.9) 


where the physical components of the strain tensor, 


e.., are given by 


ij 
0 0 
0 = 0 
Ce wales (4.10) 


For an incompressible material, the linear stress- 


Strain relation, (2.5) reduces to the ‘deviatoric relation 


= ie (Ani 
one 2u er5 ) 
where u is the shear modulus for infinitesimal strain 


from the undeformed state. 


Substituting: £On Ha from equation (4.11) in 


(4.9) we get 
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and from equation (4.10), this becomes 


Vv = »| (38) + (3) | { (4.12) 


Since e = 0 for an incompressible material, 


kk 
equation (4.10) gives the relation 

du u 

dBi a 2(¥) : (4.13) 
From equations (4.13) and (4.12) we have 


V = 6 u(2)? (4.14) 


Using equation (4.7) and the expression 


ay = Cantor : 
we find 
B 
sae = | Bama oR 
A R 


which, on integrating, gives 


2 (BoA?) 


Bo 


V = 8uTAa 
max 


Now, by equation (4.2), equations (4.8) and (4.15) 
give for the frequency, w, of the single mode of 


vibration the expression 


1 
W = Pata s + AB + i Seen) (4.16) 


Putting = N, equation (4.16) can be written as 


Ties 
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+N + 1)]? (4.17) 
and the period of oscillation, IT, is 


tT = on mB{p/[u(N? +N + 1) ] }? 


If we non-dimensionalize the period by using the 


expression 


— a | 
Lye TK De 


we find that the non-dimensional period, T, is 


bas 2 é 1, 
Tee enN (Nee --aN  +orl) canes 
For N = = = 2, the period is 2.3758. 


In the case of a thin shell, N21, and so the 
period of free radial oscillation is 1.8146 which 
agrees with the result obtained in Section 3.5 for 


the thin-walled neo-Hookean shell. 


Compressible Shell Problem in Linear Elasticity: 


The Method of Characteristics 


4.3.1 Governing equations 
The equations describing spherically symmetric 


motion are 
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J0 (o_-o,) 
EP. V: bee, ites Sts, 
PxE + yr = -2 r (4.20) 
where 
E = E/[(l+v) (1-2v)] . (4.21) 
Equations (4.18) - (4.20) form a system of three 


linear first order partial differential equations 


with O19 and v as the dependent variables and 


8 
r and t as the independent variables. These 
equations are fully hyperbolic and they permit the 
propagation of jumps in the derivatives of O.1 % 
and Vv in the characteristic directions. The 


equations of these characteristic curves can be 


derived as follows. 


Aros Eguations of the characteristic curves 


The curves in the (r,t) plane along which 
Or %% and v are continuous but their derivatives 
are not necessarily so are called characteristic 
curves and the equations governing the variations 
of O14 and v along the characteristic curves are 
called the characteristic equations. 


The characteristic curves for the set of 


@quations (4.18) - (4.20), derived in Section D.2 
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of Appendix D, Are 


+ 
Lo: os = Eo (AP22)) 


ane ey fobs Soy (4.23) 


The Ta characteristic, labelled in Fig. 4, represents 
the line along which discontinuities in the partial 
derivatives of Or Op and v propagate to the right 
with velocity c, while the i eee path of propa- 
gation of these discontinuities to the left. The 
II characteristic is the degenerate case of the 
dynamic wave at a given radial location, r = constant. 
For homogeneous materials, the quantities E, v and oe 
are constant and so the velocity of wave propagation 
is also constant. Consequently, the characteristics 
form three families of parallel straight lines. 
Substituting the expressions for dr/dt of 
equation (4.22) in equations (4.18) and (4.20) and put- 
tinge thate tor dr of equation (4.23) in (4.18)and (4.19),we 
obtain, after using the relation (4.21), the follow- 
ing characteristic equations respectively. Along te 
we have 


= v 
do. + pcdv = [-2 (0-95) ta 20C TEES, v] =, 


(4.24) 


and along II, the characteristic equation is 
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Vv 
vdo = (1-v) do, - E = dts (4.25) 


where E is Young's modulus for infinitesimal strain 


from the undeformed state. 


Aso CLONO md? SCONtInulties: in O14 T and v across 


the wavefront 


For a shell of finite thickness, say with 


A 1, the characteristic field obtained from equa- 
tions (4.22) and (4.23) is shown schematically in 
Fig. 4. The arrows indicate the path of the wave- 
front. Reflections of the wave occur on the inner 
and outer boundaries at points Cy) eb CLC mene 
is the initial propagation line and so the region ABC 
has zero stresses and velocity and is called the 
undisturbed region. 

Jumps occur in the values of Or 4 and v 
as the wave propagates back and forth in the medium 
under a step function application of internal pres- 
sure. These jumps are the difference in values of 
the quantities O14 and v ahead of and behind 
the wavefront. At the instant of reflection a double 
jump occurs, once across the Tu and the other across 
the I characteristic. 


The equations governing the propagation of 
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these strong discontinuities have been derived from 
the governing equations by CHOU and KOENIG [8] and 


these lead to the general discontinuity relations 


ey = & (4.26) 
[Oe aes oa & (4027) 
[v] = 5 = “ (4.28) 


where K is a constant which is determined from 
given boundary conditions as indicated in Section C.2 
of Appendix C. The upper and lower signs of equation 
(4.28) are for the 1 and I wavefronts respectively. 
In using the boundary condition at the inner 
surface to find K, we note that if the input function 
is not suddenly applied, K = 0, and there is no | 
discontinuity in the values of Sinan and v_ across 


the wavefronts. 


Ars. 4 Integration of the characteristic equations 


The procedure for integrating the character- 
istic equations (4.24) and (4.25), which incorporates 
the discontinuity, celations= (4626) — (42215) @1Sedis— 
cussed in Section C.l of Appendix C. The method has 
been developed from the procedure used by LEONARD and 


BUDIANSKY [ 9] to solve some beam vibration problems. 
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The response curves of the inner and outer 
surfaces of the shell with : = 2, for a step input 
pressure Q = 1.0, are shown in Fig. 19a. The curves 
indicate that this forced motion of the compressible 
thick-walled shell is non-periodic. 

The limiting case of a thin shell is examined 
by considering a suitable value of ce, the ratio of 
wall thickness to mean radius, which is given in 
terms of =, in equation (4.1), as 

(= - 1) 
Eoj= 2 Sccaape we: 
(1 + 5) 
For example, when e€ = 0.05 and Q = 0.001, we obtain 
the response curve of Fig. 19c, which shows the 
motion to be periodic. Furthermore, as Q + 0, the 


oscillation approaches a harmonic oscillation, de- 


picted by the diagram of Fig. 19d, 
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Finite Oscillation of the Incompressible Shell in 


non-Linear Elasticity 


The equation of motion of the inner surface 
of a thick incompressible spherical shell, as derived 
by ERINGEN and SUHUBI [10], is, in non-dimensional 


EOrm, 


L 


1 
ee ‘1 +(1 + 30) 3 + yee 3 - 3 ote +5) : 
x x 
4 


r SOPs - 905 
(3)p tr 5) td g (Sey), S20 


x 


(4.29) 


where, using the notation for the thick incompressible 


shell, 
ate) = ate) ; (4.30) 
5 3 
“ - (2) eile vies (4.31) 
and 
xe as 
me = dw _du 
+x 
1+y 
W r? 
where W = 7 and u(R,t) = =; .- The boundary conditions 
R 


are represented by 


Cane a te i Jie a 


Ms TR 


uK = ied fay oer | 


(Cf, >? 


Se ee et vo th il ‘iota 
PT ita cae! 


(Wes) a ome (ata - 
(16.6) of . 


they hy ' rR ? bw tye ® 


ae 
Ae RR Dae fo A sai 7 hed 
| Bi. 
| | 


as , i 


u 


Q q, (t) = q, (t) 
7 


(4.33) 


where q, (t) and q, (t) are the internal and external 


applied pressures and u is the shear modulus for 


infinitesimal deformation from the undeformed state. 


Equation (4.29) is a second order non-linear 


ordinary differential equation, the final form of 


which is obtained by substituting a given expression 


for W. 


The equation is thereafter solved numerically 


using a suitable numerical technique such as.a fourth 


order Runge-Kutta process or the method of discrete 


variables as was done for thin shell problems in 


Sections 3.2 sand 3.3. 


In the case of the neo-Hookean solid, 


(4.34) 


ee i 4 1 
dW 4,.-> ee Oe Mae oe 
du ~ 7 3 C1 u RSE gece u ) 
(47735) 
and 
3 
x —= = 
Ds, 4 () 3(1l-u —) 
yx 
1+y 


(4.36) 
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On integrating, equation (4.36) gives 


= AMG -1 -4 ei 
(say) 4x + x fe) 
aNee 
oh 
S (=) 3 (4.37) 


Substituting equation (4.37) in (4.29) we get the 


equation of motion 


(4.38) 


This last equation is solved for (i) various values 


of y, corresponding to different thicknesses of the 
shell, and a fixed value Se aconetane internal pres- 
sure Q, the results of which are shown in Fig. 23a 
and (ii) various values of Q for a fixed value of v7 
the results of which are shown in Fig. 23b. 

The response of the inner surface, or indeed 
of any given radial location, of the thick incompres- 


sible shell is found to be periodic with a period 
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and an amplitude which vary with the values of the 
internal step input pressure Q and the ratio S, of 
tremOuter=toetne snner radius. 

Pppacercurar, Lor O — 0.05, and - = 2, the 
period of oscillation of the inner surface is 2.4509. 
The response curve is shown in Fig. 23a. 

The result, plotted in Fig. 23c, for the lim- 
iting case’ ofa “thin shell for which © = 0.05 under 
a pressure Q = 0.02 is identical to the previous 
SeSuLtof Section S22 plotted in Fig. 4a, “Also, 
as Q > 0 in the thin shell limiting case, we get for 


the period of oscillation the value 1.8158 which 


agrees with the linear elasticity result quoted in 


Section? 3.2. 


Special Cases 


4,5. Thinesnelleas) limiting Casesorethe thick shell 


When € = = , where T is the wall thickness 


and R is the mean radius of the shell, is such 
that terms O(e7) are negligible, we have an approx- 


imation to the thin shell behaviour. Results obtain- 


ed from the thick shell theory then approach those of 


the thin shell analysis. For example, the period of 
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free radial oscillation for the thin incompressible 
neo-Hookean solid with c« = 0.05 is 1.8140, and the 
thick shell theory gives 1.8158. Further comparison 


of results is given in Section 6.3. 


4.5.2 Linear theory as limiting case of the 


non-linear theory of Slacercuty 


When the applied pressure Q is small enough 
to produce only vanishingly small strains, the 
constitutive equation (2.3) reduces to the linear 
isotropic stress-strain relation (2.5). Results of 


the non-linear theory are then in’ close agreement 


with the results of linear theory of elasticity. In 


particular, as Q > 0 the period of radial oscillation 


for the thin incompressible neo-Hookean shell, is 
1.8141, the value obtained from linear theory is 


1.8146. 
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CHAPTER 5 
OTHER METHODS OF SOLUTION 


For the non-linear problem of the thin in- 
compressible shell, the method of solution used by 
WANG [13] does not assume a strain energy function, 
W. Instead, the constitutive equation is written in 


the form 
= -l 
C= PI + £By+ «gB (5 <0) 


where o is the stress tensor, P an undetermined 


hydrostatic pressure, B = FFI is the left Cauchy- 


~ 


Green strain tensor and f,g are functions of the in- 


variants I I. and I. of the tensor B defined in 


1A Gxe2 3 
Secelons2 «2 

Wang derived the equation of motion under an 
internal step input pressure Q, in the notations 


used for thin shells, as 


2 2 4 4 2 
<e hOkse fo So Li eee Pet ae g 
eT da 2 4 a 2 


(5.2) 
This equation is here shown to be equivalent to 


equation (3.15) as follows. From equations (5.1) 


and (3.13), we have for Ce. and om 
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II 


9 SRE Ety EAS Ok OA (58.4) 


Using the incompressibility condition (3.14) equation 


(5.3) becomes 


Se ods ae Gace (5.5) 


Substituting the expression for o, from equation 


6 
(5-4 )eand that cor oF from equation (5.5) intom(2.12) 


we get for the equation of motion, in terms of f and 


gr 
ae8 Jota - aha) g] = (q,-q,) - ptt. 
(5.6) 

; —2 t. See a PD 7 
PMCt ing ss Ge= A] ae prt = pyEAA R and Gas ade earn Q, 
we get 

2 — 
Sep = RECO Oe eae eo 
2 3 4 
(Sica) 
which, by setting A = = becomes 
2 2 4 4 2 
iG “4 Eek ey | ee 
9-821 (5 as (= oh fe) 
Rec R i \CR i 


This is equation (5.2) obtained by Wang. 
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For the thick incompressible shell problem of 
non-linear elasticity, an alternative method, which 
does not uSe an acceleration potential, has been given 
by ZHONG-HENG and SOLECKI [11]. This method leads to 
the following form, using the same notation as in 


Section 4.4, of the equation of motion 


el 
d 5 ns ine 2 
a} 1 (2 + %) xox }, pements 255%), 


emanate) (52.9) 


which, after performing the differentiation on the 


left-hand side, becomes 


uf ee 
2x3 ee 22) = mn, Ae s-a(r +4) 


x x 
4 
= PS lL + | . ae g(x,y) = an OF 
a] eee 
C50) 


This is the dimensional form of equation (4.29). 


For the thick compressible shell problem 
of linear elasticity, a direct analytical method of 


solution is the infinite series method described by 


BAKER AND ALLEN [2]. 
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This method is based on the Midlin-Goodman technique 
which uses an auxiliary variable to eliminate time- 
dependent boundary conditions. 

The results obtained using this method are 
compared with the solutions obtained in Section 4.3 
in Figs. 19e and 19£", The graphs show that the re- 


sults are consistent. 


See ROSE, CHOU and CHOU Pele ies 
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CHAPTER 6 


DISCUSSION OF RESULTS 


Precautions in Numerical Computations 


All the programs listed at the end of this 
thesis have been written in the BASIC computer lang- 
uage for running on the 9830 Hewlett Packard desk 
CalLculaton,s, (1O,ensure Seeger cents and stability of 
the numerical computations the integration steps 
must be suitably chosen. The following integration 
steps have been found most suitable. 

(i) For the method of discrete variables and the 
fourth order Runge-Kutta process, a time in- 
terval, ot <= 001) 1s satistactorny. 

(ii) For the method of characteristics, an integ- 
ration step of 6r < 0.05 is suitable. In the 
case of a very thick shell, one for which 


eS 


A (4 
also been found satisfactory. 


10 say, a coarser step of 6r = 0.1 has 


Computations are relatively slow on the HP9830 
calculator, but the overall computing times for 
specific cases are within acceptable Limits i oO 
example, for four cycles of the response curve One 


Fig. 14a, calculated by the fourth order Runge-Kutta 
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method, using §t = 0.01, the computing time is one 
hour. For each of the curves in Fig. 19a, using the 
method of characteristics and an integration step of 
ér = 0.05, the computing time is three hours and for 
Fig. 19b, using 6r = 0.1, the total time is eight 


hours. 


Specific Results 


6.2.1 Incompressible shells 


A®periodic motion 1s obtained, for a shell 
under an internal step input pressure Q < OR 
which is the limiting value of Q discussed in 
Section B.3 of Appendix B. This type of oscil- 
lation is illustrated by phase plane diagrams and 


response curves, relating the radial stretch A to 


the time t, for the following cases. 


the The thin neo-Hookean shell for which 
€ = 0.05. Phase plane trajectories for Q = 0.05 
< Oe and Q = 0.06 > Chae y, where Ona = 0.0556 
are shown in Fig. 10a. That for the limiting 
case, Q = One t 7S shown in 10. 10D. neste 
sponse curve for Q = 0.05 is shown in Fig. 14c. 


The motion, in this case, is periodic with a 
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Period oOLeseso72. 

Free oscillation for vanishingly small 
amplitude, is represented by the phase plane 
diagram of Fig. 10d. This is for initial condi- 
tions: A, = 1.005 and A, = 0. The phase 
trajectory, in this case, is an ellipse which 
indicates a harmonic oscillation depicted by the 
response curve of Fig. 18. The period is computed 


to be 1.8140. 


The case of the thin incompressible shell 
of Mooney-Rivlin material, for which constant, 
c = 0.9 in the strain energy function of equation 
(2.6). Phase plane trajectories are shown in 
Figs. lia and llb for (1)) OF=i02060- OF ee which 
gives rise to the periodic oscillation depicted by 
Gea LS se ids)iO =) 0. 0o33e— Oe at the critical case 
andes( 111) Oe 0. 01> Oba t 

It is important to note that the periodic 
motions hitherto considered are for practically 
realisable values of stretch A. If, for instance, 
X > 3, a thin shell may rupture because it may 


not be capable of sustaining such large stretches. 


Now, an interesting case of Q > Oven is shown in 


Fig. llc where the phase trajectory is seen to 
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close again for large enough value of i, say 

A > 5. This feature is also shown for Q = ees 
in Fig. lla. Whereas these closed trajectories 
would indicate that periodic motions will occur, 


these extreme cases may not be realisable in 


practice. 


For the thin incompressible shell, of mater- 
ial with the logarithmic strain energy function of 
equation (2.7), we consider the case for which 
c= 2. Fig. 12 shows the phase plane diagrams 


fore (i) sO. = 0.0 8e< Oaat which gives rise to a 


periodic motion depicted by Fig. 16, (ii) Q 


0.0811 = Orr! the critical case and (iii) Q 
O20:3.2aa> 0. at for which the motion is non-periodic 
no matter how large a stretch i is considered. 

This last case is so because the static curve ob- 


tained for the logarithmic strain energy function 


is similar to that for the neo-Hookean form. 


The case of the thick incompressible neo- 
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Hookean shell. The response curves for forced motions 


of shells with various thicknesses, under a constant 


internal pressure Q,are shown in Fig. 23a. Those for 


a constant thickness of shell under various inter- 


nal step input pressures are shown in Fig. 23b. 
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These graphs indicate that the oscillation is 
always periodic provided Q < Ones = 0.0556. 

In the limit, as - + 1 we have the per- 
iodic response of Fig. 23c. 


Gre Zn? Compressible shells 


Hee 


The motion of the thin compressible shell, 
of Blatz-Ko material, aider an internal step fun- 
GEeLOn#application Of pressure OW< —_ is found 
to be periodic. This is illustrated by the phase 
trajectories of Figs. 13a and 13b and the response 
CUMVES FO hilgsmel/avand l/b eee inepabtaciulanaeror 
Q = 0.02, ¢€ = 0.05 and v = 0.3, the response curve 
and the phase plane diagram are shown respectively 
ine gsaelyb and sb see The peri odsomsosciVatvon;, 


fOmrCh usecase, 1S 62.69 251 


For the thick compressible shell of linear 
elasticity, the response of the inner and outer 
surfaces to an internal constant pressure are 
shown in Figs. 19a and 19b. The motion, in this 
case, is non-periodic. In particular, for - = 2 
and Q = 1, the response curves of Fig. 19a are 
obtained. For a very thick shell, for example 


one for which : = 11, the response curve of the 


, 
1h oe an he glee egeratt ris ano hom ae th 
«te? qera’ lairthan! me eed ) 

eh Le sppeeesy ae mest 
jnety sAs a Htst3fultt eh ‘eee 


iiiwos Sho bab 40i bee at) sagt Gn Gebeeteeie) oe 


107 ~pIeisniy sew wr if ina eee onet ‘“ ee ae : a 
- 


svile Refogess af (f.9 = @ tele TOP] 2 tine 9” 
Peettpgo7 Avert gt weapith anotq essdq silt has 7 
vy: tcl aden, Yo so Pieey-aet atk Be. ae weit at 

2Stea.S nt .eaetalds 708 . 


voenil to digee 
oe 

ideo te Féin) 1) So aanedian ea) , coimbhunedle : 

7 


Aye Pivepotd Ntacaate, ‘pergina ne Go cone yee : 
rT ee ae Te eee 
‘= c se8 en lyoisaad AI ciiolasenenen af) games: 

w3A W8! ..61% 3c NerEED evsoQoNTpa’) 4.4 Obes 

atenaey ami ,iioin dalds wee.o oct} §Gonienp 
ath Mol or nd wb agin A rh) ftw 1a bce 


(useo3eres AEs ait 1 ee aT 
7 


“nt 


inner surface under a pressure Q = 1 is shown in 


Fig. 19b. This shows an interesting asymptotic 
feature initiated by the shock wave at the start 
of motion. The amplitude of motion attenuates 

in such a way that the displacement approaches 
the equilibrium state asymptotically like a 
forced damped oscillation. This is because 

there is long enough cme for the imparted motion 
to die out before the shock front returns to this 
surface after reflection at the outer boundary. 
On return to the inner surface the motion starts 
again and another cycle Of oscillation is per- 
formed. 

This case of a thick shell illustrates 
what happens in the case of an infinite medium 
where the wavefront never gets reflected back 
into the medium. The response is then similar to 
one cycle of the thick shell case analogous to 
forced damped vibration. 

Another notable feature in the compres- 
sible case is the reflection of the shock front 
at the bounding surfaces of the shell. These re- 
flections appear as cusps in the response curves. 
For example, cusps occur on the response curve of 


Fig. 19a for the inner surface of a shell with 
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Sin 2 ateeaee, 46s: 

Furthermore, we note that the curves of 
Figs. 19a and 19b agree well with the response 
curves obtained by other investigators. This is 
shown by comparing the graphical results, with 
the series solution obtained by BAKER and ALLEN 


[2], schematically in Figs. 19e and 19f. 


6.3 Comparison of Results 


Results obtained by different methods are 
compared for the following cases. 
(i) Free radial oscillation of the thin neo- 
Hookean shell. 

(ii) Free and forced vibration of the thin shell 
obtained as a limiting case of the thick 
shell. 

(iii) Free and forced vibration of the thin shell 
of linear elasticity theory obtained as a 
limiting case of the non-linear theory. 

(iv) Forced vibrations of the thick compressible 
shell compared with the solutions of BAKER 


and ALLEN [2]. 
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ile Results for the period of free radial 
oscillation of the thin neo-Hookean shell, ob- 
tained by different methods are compared in 
Table 1. The answers are consistent; giving a 


periods ofsosci i lation of approximately. laclon 


Thin shell theory 


Thirds 
limiting 
case of . 
Ehwck 
shell 


C= eee nt) 
Mooney- 
Rivlin 


Direct free 


oscillation 
for Ien81-58 


a 00 Se 
he = 0 


LOVE [1] 


Analytical 


Solutions 


Rayleigh's 
eenea | 158146 


TABLE 1. Comparison of results obtained by 
different methods for the period 
of free radial oscillation of a 
thin neo-Hookean shell, for which 
€ = 0.05. 
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When € is taken small enough so that terms 

O(es) become negligible, the results of the thick 

shell theory approach those of the thin shell 

theory. Consider the following examples. 

a) Free oscillation of the neo-Hookean thin 
Shell with ¢€ = 0.05, and initial conditions: 

A, = 1.005, and h = 0. Thin shell theory 
gives a period of oscillation of 1.8140, the 
thick shell theory gives, in the es a 
value 1.8158. 

b) Forced vibration of the thin neo-Hookean shell, 
for which « = 0.05, under a constant internal 
pressure Q = 0.02. Period obtained from ore: 
thin shell theory is 2.0864, that from limiting 
Case, Orethick shel leis 2 L011 

c) Compressible thin shell, with v = 0.3, € = 
0.05. Period of free oscillation given by 
LOVE [1] is 2.3062, the result obtained from 


limiting, case of thick shelletheory 1592730547" 


As Q approaches zero the results of non- 
linear. elasticity approach those of the linear 
theory. Consider, for example, the thin neo- 


Hookean shell with ¢« = 0.05. Period obtained 


from the linear theory is 1.8146, that obtained 
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Results for the thick compressible shell 
are compared with the series solution obtained 


* 
by BAKER and ALLEN [2] in Figs. 19e and 19f. 


See ROSE, CHOU and CHOU [17]. 
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CHAPTER 7 


CONCLUDING REMARKS 


Leny DCSE O Gmiot1 On 


The non-linear elasticity problems of the thin 
and thick incompressible shells have been solved for 
admissible strain energy functions. Results indicate 
that we get the following types of response to an 
internal step input pressure, Q. Denoting by Oran! 
the maximum value of Q for which the phase trajectory 
closes, we have 

(a) LOL < Own motion is periodic with a finite 
period of oscillation, 

(ii) as Q +> 0, the response is that of linear elast- 
icity theory. Consequently, the motion is har- 
monic with a finite period of daciilatton 
Free oscillation, which involves only vanish- 
ingly small amplitudes, is a special case of 
the linear ~ theory. For this case, motion 
is also harmonic with a finite period of oscil- 
lation. 

(iii) For the linear elasticity problem of the thick 


compressible shell, where only infinitesimal 


strains are involved, the motion is non-periodic. 
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When the shell is very thick, say for ~ 20, 
each cycle of the response curve resembles 
that of forced, damped vibration. ie, ae lim- 
iting case of a thin shell, when - > 1, the 
motion approaches a periodic motion, with 
a finite period of oscillation. 

(iv) Forced and free oscillations of non-linear 
elasticity theory fou the thin compressible 


Shell are similar to those of the thin incom- 


pressible shell. 


Numerical methods 


Excellent agreement is found between the var- 
ious computed results, some closed form analytical 
solutions and the numerical results obtained by other 
investigators. This validates the results of the 
numerical methods used in this investigation. These 
numerical schemes include a modified form of the 
discrete variable method, the fourth order Runge-Kutta 
process, and the simple trapezoid rule of integration. 
Thus proven, these methods can be used reliably to 


solve spherical shell problems when analytical sol- 


utions become intractable. 
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The method of characteristics 


The method of characteristics has the advan- 
tages, over alternative numerical methods such as the 
finite-difference method, in being direct, simple and 
versatile in application. It gives good results 
especially at instants of wave relections where 
a truncated series, for example, will give only 
approximate results. Furthermore, it evaluates the 
stresses directly during the step-by-step computations, 
whereas in the alternative methods, these stresses 


have to be computed separately. 


Other input functions and boundary conditions 


Most of the results have been computed for a 
step function application of pressure but only a 
slight modification of the governing equations is 
required for other pressure and velocity inputs, for 
example, rectangular and exponential pressure inputs. 

In the calculations, non-dimensional quanti- 
ties have been used. This facilitates a direct ex- 


trapolation of the results for similar cases. 
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5. Extensions 


The following extensions can be made to the 


work presented here. 


(i) 


(eras) 


iy) 


The theories can be generalised to include 
cylindrical and plane problems by using approp- 
riate constants in the constitutive equations 
and by defining the tensor B accordingly: 

Other admissible strain energy functions, 
especially for compressible materials, can be 
investigated. Some of these, discussed by 
HADDOW and FAULKNER [12], for example the 


strain energy function 


a 5 3) - 2(vT, - 1) + erg Shy) 


[ (I OR ee 


he 
will adequately describe the response to dynamic 
loads of a useful class of compressible mater- 
ials. 

The analysis based on the method of character- 
istics can be extended to the non-linear elast- 
icity problem of the thick compressible spher- 
ical shell. This problem will be found com- 
plex because the characteristic curves are no 


longer straight lines. 
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Static pressure-stretch relation for Mooney- 
Rivlin material, material with logarithmic 
strain energy function and Blatz-Ko material. 
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Phase plane diagrams for different 
Stepeinpuc pressuuess = CUnvel (a) 
shows motion is periodic, and curve 
(b) shows motion at the limiting 
state, Of periodicity. 
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Diagram illustrating the free radial 
oscillation of a thick-walled incom- 
pressible shell: a is amplitude at 
inner radius, and u is amplitude at 
a generic radius, R. 
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Types of meshes in the characteristic 


fieLds.ot Fig. 
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Mesh» (a) 1S an 2nner 


half-mesh, mesh (b) is an interior 
full mesh, and mMesh=(¢c) 1S an Outer 


half-mesh. 
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Big. 


10a. 


Phase plane diagrams for a neo- 
Hookean thin-walled shell under 
various step function applica- 
tions of pressure. 
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Fig. 10b. Phase plance diagram illustrating 
the limiting case of periodic motion 
for a neo-Hookean thin-walled shell: 
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INPUT P,E 

PPINT TABIS"STATIC EQUILM OF BLATZ KO S.E.F." 

PPINT 

PRINT TAB25"NU="P 

PRINT 

PRINT "STRETCH","PRESSURE" 

PFINT 

S==P/ (1-2*P) 

QC=0.0C001 

PORSX=1' TO 2° STEP 0.05 

Q=00 

K=X/ (2*E) 

A=1/ (X*X) 

B=2*XK** (6*S=U) 

C=X** (12*S=6) | 

F=Q#Q/2- A* (1-K*Q) #B* ( (1K *Q) *#S) -C¥ ( (1- K*Q) 
KK {2*S= 1) ) 

F1=Q/2+A*K= B*S*K* ( (12 K¥*Q) ** (S=1)) +C* (2*S= 1) 
Pee aKa) es (25> 2)? 

Q1=Q0-F/F1 

IF ABS(Q1-Q0) <= 1E-10 THEN 85 

QC=Q1 

GOTO 30 

PRINT X,Q1 

Qc=01 

NEXT X 

END 


Poogran el. serogram tor solving the equilibrium 
equation for the thin compressible 
shell of Blatz-Ko material 
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C,E,B 
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ee 
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Scete | WREME fhe *y 


PF Gt98¢2<0 THEN $5 
7 if TeTI=T THEY 260 
TF % >= & THEW DT 
Gi=G) 
g=o+n 


AS GOSUP 25° 
320 Pier 


; E Be .0te £ ar on fave tof metpoxrt 

a3 Sth SF loddbg TO sy! sSmeastq Jugeti 

J-HSNGOM (TS siantiinicns noes ts to 
f=-X3 . -isiteatsm alivis 

Geson 23¢ / 

F3«F 

7 ee SBS (kteR) <= 18190 TREW 210 

=. 2) Pier2<¢ inru 80 

© fF Fier3>° Tere 195 

TP PieeseO tery 210 

© 2253 

283 Wk? 

E#90 Goto 125 

a5 s4=x5 

200 Rex! 

GOTO 1796 

‘21C PPINT 

215 Print 


"OQ (RA) FOF DTP. OTAB.CY A GT¥EN 5.8.9." 


-S£ msxpord 


 SPANT TABS CEITICeL STSeTCH=""43,"%O (MAX) =" 


30 82 


Be teed) 72 17 (a (8848) 
oom dail lai 


she (3043-1) /> SeC*#te D4 (VC) Oe 372 


60 P 


; 0 
oi EP $06 tare 225 


f 
ome 
a 


pedece (4/t~ 1/ (Re*T)) +2" (1-C) * 17 (OSD) 


bef "40 OTKEN POUT UP £¢/k UP TO x=*x 
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oa se | mY 


Program 2a. Program for evaluating the maximum step 


input pressure, Os for periodic motion 


of an incompressible shell of Mooney- 
Rivlin material. 


INPUT C,E,H 


PRINT "Q(MAX) FOR DYN.STAB.OF A GIVEN S.E.F." 


PPINT 

PRINT HE="C ; "E="E, WH="H 
PRINT 

G1=0 

X=1+H 

GOSUB 230 

G2=F 

FESGI*GZ< 0 THEN 95 
IF INTX=X THEN 260 
pioeke2= oo THEN) 276 
G1=G2 

X=X+H 

GCTO 55 

X1=X-H 

X2=X 

R=X2 

X=X1 

GOSUB 230 

F1=F 

X=X2 


! GOSUB 230 


F2=F 

X3=X1= (X2=X1) *F1/ (F2-F 1) 

X=X3 

GOSUB 23C 

F3=F | 
IF ABS(X3-R) <= 1E=10 THEN 210 
IF F1*F3<0 THEN 180 

IF F1*F3>0 THEN 195 

IF F1*F3=C0 THEN 210 

X2=X3 

R=X2 

GOTO 125 

X1=X3 

R=X1 

GETON 110 

PRINT 

PRINT 


PRINT TABS"CRITICAL STRETCH="X3,"Q (MAX) ="Q 


END 


A=2*C* (1/X-1/ (X**7) ) +2* (19C) * (X= 17 (X**5) ) 


B= (X**2) /2+1/ (4* (X**4) ) 
D= (X**4) /4+1/ (2* (X¥**2) ) 
Q=E*A 


F=A* (X*%* 3-1) /3> 2*C*B= 2* (1°C) *D+3/2 


RETURN 


PRINT "NO OTHER ROOT OF X'/X UP TO X="X 


GOTO 80 

IF H<1E-06 THEN 225 
H=H/10 

GOTO 35 
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Program 2b. Program for evaluating the maximum step 


input pressure, Daas? for periodic motion 


of an incompressible shell of material 
with logarithmic strain energy function. 


VA de) 
270 
210 
280 


INPUT C,E,H 
PRINT "Q(MAX) FOR DYN.STAB.OF A GIVEN S.E.F." 
PPINT 

PRINT "C="C,"F="E, "y= 0H 

PRINT 

G1=0 

X=1+H 

GOSUB 230 

G2=F 

IF G1¥*G2<0 THEN 95 

IF INTX=X THEN 260 

DRox) 2= 5 THEN 270 

G1=G2 

X=X+H 

GOTO 55 

X1=X-H 

X2=X 

R=X2 

X=X1 

GOSUB 230 

F1=F 

X=X2 

GOSUB 230 

F2=F 

X3=X = (X2=X1) *F1/ (F2°F 1) 

X=X3 

GOSUB 230 

F3=F 

IF ABS(X3-R) <= 1E-10 THEN 210 

IF F1*F3<0 THEN 180 

IF ¥1*F3>0 THEN 195 

IF F1*F3=0 THEN 210 

X2=X3 

R=X2 

GOTO 125 

X1=X3 

R=X1 

GOTO 110 

PRINT 

PRINT 

PRINT TABS"CRITICAL STRETCH="X3,"Q (MAX) ="Q 
END 

A=2¥*C¥ (1/X=1/ (X¥*7) ) + 3% (1-C) * (X**3= 1/7 (X** 3) ) 7 (24 X¥**6) 
B= (X**2) /24+1/ (4% (X¥**4) ) = 3/4 

D= (LOG ( (2+X**6) /3) ) /6= (LOG ( (3*X**6) / (2+ X**6) )) /12 
Q=E*A 

F=A* (X**3= 1) /3> 2*C*B= 3* (1°C) *D 

RETURN 

PRINT "NO OTHER ROOT OF X'/X UP TO X="xX 
GOTO 80 

IF H<1E-06 THEN 225 

H=H/5 

GOTO 35 
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INPUT Q,E,P,H 
PRINT BO="O, NE=tE, "NO="p, WH = tH 
BEINT 

sabe "Q(MAX) FOR BLATZ KO MATERIAL" 
je 

x0=1 

VC=0 

DO=1 

C==-1/ (1-2*P) 

GOSUB 50 

X=X0 

GOSUB 45C 

xXC=X 

N=C 

AO= (Q*¥XO*X9/E) -2* (X0-D0**2/X0) 
A1=A0 
V1=VO+(A1+#A0) *H/2 
X1=X0+(V14#V0) *H/2 

X=X1 

GOSUB 450 

A2= (Q¥X*X/E) = 2* (X=D** 2/X) 

Pigea hs (AZ@A 1) Cle=O8 THEN» 11/5 
A1=A2 

GOTOR/5 

DRA AZ<0 THEN 500 

N=N+1 

vc=Vv1 

XC=X 

DO=D 

IF N*H >= 500 THEN 750 

GOTO 55 

RETURN 

D=D0 

F=Q4+ (2*D/ (X*X)) = (2* (X¥*X#D) *¥C) 
F1= (2/ (X*X) ) = (2*C* (X** (2*C) ) ¥(D** (C-1))) 
D1=D0= (F/F1) 

IF ABS(D1-D0) <= 1E=10 THEN 485 
DO=D1 

GOTO 450 

D=D1 

RETURN 

IF J=2 THEN 600 

J=J3 +1 

GOTO 120 

L=1 

Xx0=X 

vo=vV1 

DO=D 

IF L=3 THEN 650 


> H=H/2 


L=L+1 

GOSUB 50 

IF VO <= 1B=08 THEN 700 
IF vO>0 THEN 680 

PRINT "APPROX Q (MAX) ="0 


Dee POR) ¢ [DVL O*8) 82°C) (4 


a) 
a. wwe ‘ , 
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665 
670 
680 
685 
690 
700 
750 


Q=O+ (0/10) 
GOTO 35 


PRINT "APPROX Q (MAX) ="¢Q 


Q=Q=- (0/10) 
GCTO 35 


PRINT "CRITICAL STRETCH="X0,"Q (MAX) ="0 


END 


Program 2c. 


Program for evaluating the maximum 


step input pressure, Oa aah for per- 


iodic motion of a compressible shell 
of Blatz-Ko material 
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FORMAT MIX Eo OX 229), 6X ph 12. 96x aE 2.9 
FORMAT MX, SOR, 4k 76 6.3,4X, F657 0x 9 ou 
ENEUTISO,&, C,H 

PRINT 

WRITE (15,12) "Q="Q,"E="E, "C="C, "y=" 

PFINT 

PRINT TAB11"DYNAMICS OF A COMPR THIN SPH SHELL" 


PRINT 

PRINT TAB11"METHOD 2: MODIFIED DISCRETE VARIABLE" 

PRINT 

PRINT TAB11"TIME"TAB24"ACCN" TABU I"STRETCH "TAB5S6 
"RATE OF STRETCH" 

PRINT 

PRINT TAB12"0.0"TAB4O"1. O"TABS8"0.0" 

x0=1 

VC=0 

O=1 

N=0 

X=X0 

GOSUB 450 

AC=A 

A1=A0 

V1=V0+(A1+A0) *H/2 

X1=X0+ (V14V0) *H/2 

X=X1 

GOSUB 450 

A2=A 

IF ABS(A2=-A1) <1E=08 THEN 120 

A1=A2 

GOTO 75 

N=N+1 

VO=V1 

X0=X 

T=N*H*O 

WRITE (15, 10) T,A2,X0,VO 

IF T >= 500 THEN 180 

O=0+1 

GOTO 55 

END 

R=Q¥X*X/E= 2*C* (X= (1/X¥**5) ) = 2* (TC) * (X#*3= (1/7 X** 3) ) 

RETURN 


Program 3a. Discrete variable method for solving 
the governing equation for a thin in- 
compressible shell of Mooney-Rivlin 
material 


ee ete ee’ 


aie ADOT IPP RST 


Haauey, 


ee tires Sor 444 
el & fo nei ys 
| ag sasay re 


2 108 
a 7 . . 


. TIE, PO. 2,717. 9,6%, 99.9, 68, 818.9 
} 2 V18, 86, 2,03, 76.3,48%, 76,9, 6077.8 
’ 5 Q,%,¥,8 


© PE 
7 ITE €15, 42) "O=%q, "E="£, "Cie ®p, het, 


* PRINT TAOVITUYHARICS OF A CoapT rete aon Sein = 
‘24 PEINY 
ee SELST TABI“ HETHOD 2: KODIFIED orseeerE Viet ASiEe 
; PRINT 
PESST TAD TIME*TAD24*ACCH*TARG et RECN TARSE 
*, "RATE Of STRETCH" 

32 PRINT 
oe PRINT TAL IDO, Orranahes. ower, 
40 XC=<1 

92 WH=0 


a .485 Ob=' 


SS Cee t/ (\- 24 2) 
>» S50 Xero 
: SZ GOS0B 45° 
|, 83 x0=2 


$4 O<1 
, $55 t-0 
ov set tae Ostrom 41Usse8e7ANGrosio df métpord 
79 BAZAY @ 103 nOs¢eypsS onimzsvop ars 
7S PieVOs4+hs ehG4O/ tise of dia2zsxamos 
80 BISAGe (V1+' ous; [sizetem 
».85 yer) } 
GOSOR 4>G6 
Ag@ t77kOt/ Fr) -2* (i=pe*cst 
ZY AHS (A2D~? 
95 £3=42 
70 GoTo 75 
920. Reyes : 
v= 7 
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DCepD 7 
Tf #eteO, 1 Tihew 13% 
GOTO 65 
T=8*h +O 
SETTER (75, 70) 7,482,480, 7 
1? T >= 10 THEN TAC 
O=pH+1 
GOTO 55 
> Beep 
O50 DebOd 
we Pee e2eD/s Feri) © (2° (TFS FD) FF) 
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Program 3b. Discrete variable method for solving 
the governing equation for a thin 
compressible shell of Blatz-Ko 
material. 
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FORMAT 11X,F4.1,6X,F12.9,6X,F12.9,6X,F12.9 
FORMAT 11X,F6.3,4X,F6.3,4X,F6.3,4X,F7.4 
INPUT Q,E,P,H 
PRINT 
WEITE (15,12) "Q="0,"E="E, "NU="P, "H="H 
PRINT 
PRINT TAB11"DYNAMICS OF A COMPR THIN SPH SHELL" 
PRINT 
PFINT TAB11"METHOD 2: MODIFIED DISCRETE VARIABLE" 
PRINT 
PRINT TABI1"TIME"TAB24"ACCN"TABGI"STRETCH"TABS6 
"RATE OF STRETCH" 
PRINT 
PRINT TAB12"0. O"TABGO"1, O"TABS8"0. 0" 
X= 
vVo0O=0 
DO=1 
==1/ (1=2*P) 
t=x0 
GOSUB 450 
XG=x 
O=41 
N=0 
AO= (O*X0*XO/E) = 2* (X0=D0**2/X0) 
A1=A0 
V1=V0+(A1+A0) *H/2 
X1=X0+ (V1+V0) *H/2 
Y=X1 
GOSUB 450 
A2= (Q*X*X/E) -2* (X= D**2/X) 
IF ABS (A2=-A1) <1E°08 THEN 120 
=A 
GOTO 75 
N=N+1 
vVC=V1 
XC=X 
DO=D 
IF N¥*H=0.1 THEN 130 
GOTO 65 
T=N*H*O 
WRITE (15,10) T,A2,X0,V0 
thet >= 10 THEN 180 
O=0+1 
GOTO 55 
END 
D=D0 
F=Q+ (2*D/ (X¥X)) > (2% (X¥*X*D) F*C) 
F1=(2/ (X*X) ) = (2*C* (X** (2*C) ) * (D** (C=1))) 
D1=D0= (F/F1) 
IF ABS(D1=D0) <= 1B=10 THEN 485 
DO=D1 
GOTO 450 
D=D1 
RETURN 
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Program 4a. Fourth order Runge-Kutta process for 


solving the governing equation for a 
thin incompressible shell of Mooney- 
Rivlin material. 


FORMAT 11X,F4.1,6X,F12.9,6X,F12.9,6X,F12.9 

INPUT Q,E,C,H 

PRINT 

PRINT "O="g, WR= An ALE Ol LT ,"H="9 

PRINT 

PRINT TAB11"DYNAMICS OF A COMPR THIN SPH SHELL" 

PRINT 

PRINT TAB11"METHOD 1: 4TH ORDER RUNGE KUTTA" 

PRINT 

PRINT 

PFINT TABI1"TIME"TAB24"ACCN"TABUI"STRETCH "TABS6 

"RATE OF STRETCH" 

PRINT 

PRINT TABI2"0.0"TAB4O"71.O0"TABS8"0.0" 

XxC=1 

Vo=0 

DO=1 

O=1 

N=0 

X=X0 

GOSUB 500 

K1=K 

X=X0+H*V0/2+H*K 1/8 

GOSUB 500 

K2=K 

X=XO#+H*VO+H*K2/2 

GCSUB 500 

K3=K 

X1=XO+H*VO+H* (K1+2*K2) /6 

V1=V0+ (K1#4*K2+K3) /6 

A1=Q*X1*X1/E= 2*C* (X1= (1/X 1**5) ) —2* (1-C) * (X1**3= 
(1/X 1**3) ) 

N=N+1 

XC=X1 

Vo=V1 

IF N*¥H=C0.1 THEN 140 

GOTO 60 

T=N*H*O 

WRUTE (15,10) T,A1,X1,V1 

IF T >= 10 THEN 180 

0=0+1 

GOTO 55 4 

END 

K=H* (Q®X¥X/E= 2*C* (X= (1/X**5) ) —2* (19C) * (X*#3- 
(1/X**3))) 

RETURN 
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FORMATS 11X,F01,6X,F 12.9,6X,F12.9,6X,F 12.9 

FORMAT 11X,F6.3,4X,F6.3,4X,F6.3,4X,F7.4 

INPUT Q,E,P,H 

PRINT 

WRITE Gis) NO="Q0, "E="E,"NU="P, tH="H 

PRINT 

PRINT TAB11"DYNAMIC® OF A COMPR THIN SPH SHELL" 

PPINT 

PRINT TAB11"METHOD 4: 4TH ORDER RUNGE KOTTA" 

PRINT 

PRINT 

PRINT TAB11"TIME"TAB24"ACCN"TABGI"STRETCH"TAB56 
"RATE OF STRETCH" 

PRINT 

PRINT TAB12"0.0"TABGO"1, O"TAB58"0.0" 

XO=it 

vo=0 

DO=1 

C==1/ (1=2*P) 

O=1 

N=0 

X=xX0 

GOSUB 450 

GOSUB 500 

K1=K 

X=XO+H*VO/2+H*K1/8 

GOSUB 450 

GOSUB 500 

K2=K 

X=XO+H*VO+H*K2/2 

GOSUB 450 

GOSUB 500 

K3=K 

X1=XO+H*VO+H* (K1+2*K2) /6 

V1=VO+t (K1+4*K2+K 3) /6 

A1= (Q*X1*X1/E) = (2*X1) + (2*D*D/X1) 

N=N+1 

X0=X1 

vc=Vv1 

DO=D 

IF N¥H=0.1 THEN 140 

GOTO 60 

T=N*H*O 

WELTE) (15,10) Dy Ad, x lev) 

TheTe>=910) THEN) 180 

O=0+1 

GOTOeD> 

END 

Bene (2* (X*X*D) **C) 

F=Q+ (2*D/ (X*X))- 

ae oy reek) y= (2*C* (X¥* (2*C) ) * (D** (C1) ) ) 

DI=D08Fr/F I 

IF ABS(D1-D0) <= 1E~10 THEN 485 
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475 DC=D1 

480 GOTO 450 

485 D=D1 

490 RETORN 

500 K=H* (Q¥X*X/E= 2*X+2¥*D*D/X) 
510 RETURN 


Program 4b. Fourth order Runge-Kutta process for 
solving the governing equations for 
a thin compressible shell of Blatz-Ko 
material. 
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INPUT Q,E,C,H 
T=0 

X0=1 

vC=0 

N=1 

X=X0 

L1=V0 

GOSUB 400 

AOQ=A 

A1=A0 

V1=V0+(A1+#A0) *H/2 

X1=X0+ (V14+V0) *H/2 

X=X1 

GOSUB 400 

A2=A 

IF ABS(A2=A1) <1E-08 THEN 40 
A1=A2 

GOTO 24 

IF N/10=INT(N/10) THEN 80 
IF L1*V1<0 THEN 54 

N=N+1 

VO=V1 

XO=X 

GOTO 14 

PRINT 

T=T+ (N=1) *H 

PRINT TAB1S“APPROX PERIOD OF OSCILLATION="" 2*T 
PRINT 


) IF H=1E=10 THEN 74 


H=H/10 

VO=L1 

XC=X 

GOTO 12 

END 

PRINT T+N*H,V1,X 

GOTO 42 

A=Q¥X*X/E= 2*C* (X= 1/ (X**5) ) = 2* (V9C) * (X**3= 17 (X**3) ) 
RETURN 


Program 5a. Program for finding the period of 
oscillation of a thin incompressible 
shell of Mooney-Rivlin material 
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INPUT Q,E,C,H 
T=0 


L1=v9 

GOSUB 400 

AOQ=A 

A1=A0 

V1=V0+(A14A0) *H/2 

X1=X0+ (V14V0) *H/2 

X=X1 

GCSUB 400 

A2=A 

IF ABS(A2-A1)<1E=08 THEN 40 
A1=A2 

GOTO 24 

IF N/10=INT(N/10) THEN 80 
IF L1*V1<0 THEN 54 

N=N+1 

VC=V1 

XO=X 

GOTO 14 

PRINT 

T=T+ (N= 1) *H 

PRINT TAB15"APPROX PERIOD OF OSCILLATION="2*T 
PRINT 

IF H=1E-°10 THEN 74 
H=H/10 

VO=11 

XO=X 

GOTO 12 

END 

PRINT T+N*H,V1,X 

GOTO 42 

B1=X= (1/ (X**5) ) 

B2=X**5= (1/X) 

B3=2+ (X**6) 
A=Q*X*X/E= 2*C*B 1= 3* (1=C) *B2/B3 
RETURN 


Program 5b. Program for finding the period of 


oscillation of a thin incompressible 
shell of material with logarithmic 
strain energy function 
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Program 5c. Program for finding the period of 
oscillation of a thin compressible 
shell of Blatz-Ko material. 


DC=1 
INPUT Q,E,P,H 

C==1/ (1-2*P) 

T=0 

X0=1 

vV0=0 

N=1 

X=X0 

L1=Vv0 

GOSUB 400 

AC=A 

A1=A0 

V1=VO0+t (A1+A0) *H/2 
X1=X0+4(V1+V0) *H/2 

X=X1 

GOSUB 400 

A2=A 

IF ABS(A2=A1)<1E-08 THEN 40 
A1=A2 

GOTO 24 

IF N/10=INT(N/10) THEN 80 
IF L1*V1<0 THEN 54 

N=N+1 

VO=V1 

XO=X 

GOTO 14 

PRINT 

T=T+ (N=1) *H 

PRINT TAB1S"APPROX PERIOD OF OSCILLATION="2*T 
PRINT 

IF H=1E-10 THEN 74 

H=H/10 

VO=L1 

XO=X 

GOTO 12 

END 

PRINT T+N*H,V1,X 

GOTO 42 

GOSUB 450 
A=Q¥X¥*X/E=2* (X= D**2/X) 
RETURN 

D=DO 

F=Q+ (2*D/ (X*X) ) = (2% (X*¥X*D) **C) 
F1=(2/(X*X) ) = (2*C* (X#* (2*C) ) *(D#* (Co 1))) 
D1=D0= (F/F1) 

IF ABS(D1=D0) <= 1£-10 THEN 485 
DO=D1 

GOTO 450 

D=D1 

RETURN 
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PRINT "Q=APPLIED PRESSURE","P=POTSSONS RATIO" 

PRINT "Z=SHELL THICKNESS","K=INTERVAL SIZE" 

PRINT "M=NUMBER OF VERT.GRID LINES=(Z/K) +1" 

PRINT "R9=PLOTTING RADIUS", "T6=INITIAL TIME" 

PRINT "INPUT Q,P,Z,K,M "“,"INPUT R9,T6" 

INPUT Q,P,Z,K,M 

INPUT R9,T6 

C1=2* (1=P) /(1-2*P) 

C2=P/(1-P) 

FORMAT F7.3,2X,F6.3,2X,E15.8,4X,£15.8,4X,E15.8 

DIM H(105) ,U(105) ,#(105) ,P(105) ,Q (105) ,R (105) 

DIM A(3,3),S(3,1) ,B(3,1) 

DIMC(2, 2) ,D(2, 1) ,E (2,1) 

Dimer (2,2) 2 (2, 1)76 2.4) 

PRINT "CHARACTERISTICS METHOD (TRIANGULAR REGIONS) " 

PRINT 

PRINT "APPLIED PRESSURE="0,"POISSON'S RATIO="P 

PRINT "SHELL THICKNESS="Z,"INTERVAL SIZE="K 

PRINT 

READ RO,TO 

DATA 1,0 

PRINT TAB1"TIME"TABO"RADIUS"TAB21"S(R) "TABU2"S(T) " 
TABS7"PARTICLE VEL" 

GOSUB 952 

T=1 

I3=14=1 

M1=M 

A=1 

IF INT(I/2)=I/2 THEN 290 

R=R0 


1 FOR N=1 TO M 


IF I#1 THEN 130 

P(N) =Q (N) =R (N) =0 

H (N) =P (N) - (I3*Q/R) 

U (N) =O (N) = (13*C2*Q/R) 


) W(N) =R(N) + (13*Q/ (C1*R) ) 


R=R+K 
NEXT N 

I3=-1%13 

GOSUB 920 
T2=TO+ ( (I= 1) *2Z) 
R2=RO0 

I1=1 

N=1 

GOSUB 630 

GOSUB 685 

IF I11=Z/K THEN 280 
N=2 

T1=T4 

R1=R4 

GOSUB 770 

GOSUB 410 

IF N=M1 THEN 250 
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N=N+1 

GOTO 215 

T2=T4= ((M1=1) *K) 
R2=RO 


) GOSUB, 920 


T1=11+1 

GOTO 185 

GOSUB 920 

IF I=39 THEN 285 
I=I+1 

GOTO 95 

PEN 

END 

R=RO+Z 

FOR N=1 TO M 

H (N) =P (N) + (T4*Q/R) 
U(N) =Q(N) + (14*C2*Q/R) 
W (N) =R(N) + (I4*Q/ (C1*R) ) 
R=R-K 

NEXT N 

T4s-14I4 

GOSUB 920 
T2=TO0+((I=1) *Z) 
R2=RO+Z 

I2=1 

N=1 

GOSUB 1005 

GOSUB 545 

IF I2=Z/K THEN 280 


R1=R4 

GOSUB 1065 

GOSUB 410 

IF N=M1 THEN 375 
N=N+1 

GOTO 345 
T2=TU= ( (M1-1) *K) 
R2=RO0+Z 

GOSUB 920 
I2=12+1 

GOTOs315 

Z3=R4/K 

Z1=R1/K 

Z4=R3/K 

A (1,1) =1+23 
A(1,2)=<-1 

A (1,3) =C1* (“C2=-2Z3) 
A (2,1) =2Z3-1 
A(2,2)=1 

A (2,3) =C1* (23-C2) 
A(3,1)==-P 
BS 72) six 2 


VTS 


A (3,3) =~2* (14+ P) /Z3 

IF INT(I/2)=I/2 THEN 480 
B(1,1)=X1* (Z1=1) #¥1#V1*C1* (C2=71) 
B (2, 1) =X3* (1+ 24) -¥3+V3*C1* (Z4+C2) 
GOTO 490 

B (1,1) =X3* (Z4=1) +¥3+V3*C1* (C2=-Z4) 
B (2,1) =X1* (1421) = Y¥14V1#C 1% (Z14C 2) 
B(3, 1)=-X2*P+Y2* (1-P) #V2*2* (14+P) /Z3 
MAT A=INV(A) 

MAT S=A*B 

X4=S (1,1) 

Y4=S (2,1) 

V4=S (3,1) 

H (N) =X4 

U(N) =Y4 

W(N) =V4 

RETURN 

X4=0 

Z3=R3/K 

Z4=R4/K 

F(1,1)==-1 

F (1,2) =C1* (-C2=Z4) 

F(2,1)=1=P 

F (2,2) =-2* (1+P) /Z4 


G (1, 1) =X4* (< 1-24) #X3* (Z3=1) #¥3+V3*C1* (C2= 23) 
G (2, 1) =X4*P=X 2*P+Y2* (1=P) #V2*2* (1+P) /Z4 


MAT F=INV(F) 
MAT T=F*G 
Y4=T (1,1) 
V4=T (2, 1) 
H (N) =X4 
U(N) =Y¥4 
W(N)=V4 
RETURN 
X2=H (1) 
Y2=U (1) 
V2=W (1) 
X3=H (2) 


) ¥3=U (2) 


V3=W (2) 

T3=T2+K 

R3=R2+K 
TH=T2+2*K 

R4=R2 

RFTURN 

X4==0 

Z3=R3/K 

Z2=R4/K 

orb 

@(1) 2) =C1* (Z2-C2) 
C(2,1)=1-2P 

C (2,2) =-2* (1+P) /Z2 


E (1,1) =X4* (1-22) +X3* (14+23)-¥3+V3* (Z3+C2) *C1 
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E (2,1) =X4*P=X2*P+Y2* (1=P) +V2*2* (14P) /Z2 
MAT C=INV(C) 

MAT D=C*E 

Y4=D (1,1) 

V4=D (2,1) 

GOSUB 962 


H(N) =x4 
U (N) =¥4 
W(N) =V4 
RETURN 
J=N=1 

L=N+1 

X1=H (J) 
Y1=0 (J) 
V1=W (J) 
X2=H (N) 
Y2=U (N) 
V2=W(N) 
X3=H (L) 
Y3=U (L) 
V3=W (L) 
T2=T1-K 
R2=R1+K 
T3=1T1 


R3=R1+2*K 


TOU=T1+K 
R4=R2 
RETURN 
N=M1 


P (A) =H (N) 
Q (A) =U (N) 
R (A) =W(N) 


A=A+1 
M1=M1=1 
RETURN 


SCARUELO,40,-0.2,0. 7 
XAXIS 0,2 
WAR RECS ee Ulg O'en! 


FO=0 


PEROT. 16, 0 


RETURN 


F1=R4* (Y4* (1=P) -P¥*X4) / (2* (1+P) ) 
WRITE (15,20) T4,R4,F1,Y4,V4 
PLOT T4,F1 


RETURN 
T3=T2+K 
R3=R2-K 


TU=T2+ (2*K) 


RU=R2 
L=N+1 
X2=H (N) 
Y2=0 (N) 
V2=W (N) 
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Program 6. 


X3=H (L) 
Yo-U(L) 
V3=W (L) 
RETURN 
J=N=1 

L=N+1 

X1=H (J) 
Y1=U (J) 
V1=W (J) 
X¥2=H (N) 
Y2=U (N) 
V2=W (N) 
X3=H (L) 
Y 3=0 (L) 
V3=W (L) 
T2=T1-K 
R2=R1-K 
T3=71 


R3=R1= (2*K) 


TU=T1+K 
R4=R1-K 
RETURN 


Method of characterists for solving 
the governing equations for the 
thick compressible shell 


22 


5 PRINT “"NEO-HOOKEAN SPHERE STEP FUNCTION APPLICATION 
OF PRESSURE" 
10 DISP "Q,N"; 
20 INPUT Q,N1 
SCSULSDBUINT (STEP Si7 reat. 
40 INPUT T 
50 PRINT "Q="Q,"N="N1,"INT STEP SIZE="T 
COMmSCALESO, So 1014 
7TOMXAXTS 1,025 
®0 YAXIS 0,0.925 
90 N=0 
100 G=N1**3=1 
Pie Yess) 
120 V=0 
130 PLOT N¥*T,X 
135 IF N¥T/O. 1#INT(N*T/0.1) THEN 150 
140 PRINT "T="N*T,"X="X,N yay 
150 D=(1+G*X** (=3) ) ** (=1/3) 
160 E=((GtX**3) /(1+#G) ) ** (= 1/3) 
170 F=G*X** (= 3) 
18C A== (V*¥V* (3= (F*¥D**4) /(1=D) ) = (4/X4X** (24) =4¥E= E**U 42 *0) 
4 (1-D))7 (2*5) 
190 A2=A 
200 W=V+(A2+A) *T/2 
210 Y=X+ (W+V) *T/2 
ee CU= (SL EG*Y S¥ (2) a) ei(= 1773) 
230 E=((G+Y**3)/ (1+G) ) ** (= 173) 
240 F=G*Y** (=3) 
250 Al== (WEW* (3= (F¥D*¥*4) / (12D) ) = (4/4 Y** (<4) = 4*E= EX*4 +2*Q) 
7 AtSD):) 7 (2528 
260 IF ABS(A1~A2) <1E=08 THEN 290 
276 22=h4 
280 GOTO 200 
290 N=N+1 
310 X=Y 
320 V=W 
330 IF N*T=5 THEN 350 
340 GOTO 130 
350 PEN 
360 END 


Program 7. Program for solving the governing 
equation for the thick-walled 
incompressible neo-Hookean shell 
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APPENDIX A 
NUMERICAL METHODS AND ERROR ANALYSES 


Numerical Scheme for the Method of Discrete 


Variables 


For a time interval At > 0, let ty =k At, 
k = 0,1, +--+ be the time after the kth interval, and 
let a particle in motion be at position hie at time 
ty: Then, we een define the particle velocity 
i = 1 (Cele ke= 12, aa cand,thesparticle SS ee 


lon hye =a (ty), k = 1,2, ... at time t, by the 


equations 
i] | 
he + hyn] _ he - yey hae 
2 At v et 
and Gert 
1] w q 4 
A, +d =e 
k k-lg an ok k-1 Shea) Ona 
2 At 
(A.2) 
respectively. 


If the particle acceleration, Mies SWE} Gh vaste te 
ction 
ier (3) 


x a N Oy 


then, an improved version of discrete mechanics can 


be described by the following algorithm, which is 


eG 
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peels menage, wo 2 . 
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written in BASIC computer language in Program 3a. 


Ci eeeSGot ke 1n 


it 
(ii) Read Mend and Mend 


(SUaIaLy! Gs aie! Ayo from equation (A.3) . 


(iv) w x w 
Set he eed > 


(v) Find an approximate Me using equation (A.2), 


and from this, an approximate A using 


pou 
equation (A.1). 


(vi) Substitute the approximate \, in equation 


k 


k 
(A.3) to find an approximate A 


Ly oe 


, denoted by 


(Vil) Transfer A ianto hie 3 


(viii) Repeat steps (v) to (vii) until hy is found 
within a specified error tolerance, e. 


(ix) Substitute the final A, into equation (A.2) 


k 


, and this in equation 


k 
(A.1) to find the final A 


to glind thesia nalen 
er 
' ' 
(x) Transfer A, into Ay,z and », into Ay-1 0 
proceed to the next time interval, tea] A 
(xi) Repeat steps (iii) to (x) for the interval, 


esl 


The flow chart for this algorithm is shown 
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Error Analysis of the Discrete Variable Method 


Equations (A.1) and (A.2) are basically the 
trapezoid rule of integration for which the error, 


E, is given by 


h? w 
Er - TD fe Aa) (A,4) 
where h = At is the constant time interval which 


would be specified, and n is such that 0 < n < h. 


Then, for the particle velocity, v the error 


ved 
3) 
a5 a(n) (A.5) 


which by forward finite difference formula of dif- 


ferentiation becomes 


Sageiec die) sakaeke Loy (te ea a 
iy h 12 k k-l 
(A.6) 


where (a,-a,_4) is the maximum difference, over the 
whole time range, between the particle accelerations 
at consecutive times t,_1 and th. 

Simi larily, =the error E. in evaluating the 


particle displacement Xr is 
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where (v,- V,.,) is defined in a similar way to 


a ag i 


The values of (a,- a ) and is ) can 


viet 


be found in specific cases. For example, if we take 


Kel 


h = 0.01, then, for a thin neo-Hookean shell with 

Ee = 2 = 0.05, under an internal step input pressure 
OF=90..0:5, (ay- ay_3) = 0.021 and (as Ae = 0.4849. 
Hence, from equations (A.6) and Cay Ej = — 0/175 x 
7 Ome and LA ene 0.4041 x Wi The specified error 
tolerance, E, on the acceleration a(t,) is 

- 0.01 x Tome Consequently, the total error, Ey! ral 


the computed particle displacement x(ty), is 


t = 
Pe oor Ger re ot o8seol Sloe: 


(Oye 
ba ortogie ee antes 


: 3 
In general, since Ey and E, are O(h )\,, .and 


: 3 
E can be chosen within Otho then E,»is also sGtn=). 


The Fourth Order Runge-Kutta Process 
Let the acceleration  , be a function 


X =A, (A) (A.8) 
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of the stretch A, then the fourth order Runge-Kutta 
process uses the following iterative formulae to 


solve the differential equation (A.8). 


= ; Low? 
atl Ns + hho + ra h (K, + K. + K,) 
(A.9) 
t i] a 
bein = ee + gn (Ky + 2K, + 2K, + K,) 
(A.10) 
where 
Ky =A (A_) 
hs W h t 
K, =r (A. + 5 aes) 
Be 
ae byt yb 
ee ne ee 
WW t h2 
oy = r (A, a‘ ne ge Dic Ky) (Aer) 


and h = At is the constant time interval between pos- 


itions represented by displacements dn and ntl : 


1 
and Ant he are respectively the displacement and 


velocity at time to: A program using this method 

to integrate the differential equation governing the 
motion of a thin incompressible shell is shown in 
Program 4a. 


The error, E, in an ret order Runge-Kutta 


method is given by 
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(2) (1) 
ye = 
E . A (Arie) 
26%= 11 
1 
where va ) and , (2) are the values of A at time th 


using step sizes h and 2 respectively in stepping 
forward from t.-1: Values for these can be found in 
specific cases. For example, in the forced vibration 
of the thin neo-Hookean shell for which e€ = 0.05, 
under an internal step input pressure Q = 0.055, 


1) = 1.561290868 ana 1'7) = 1.561289352. Since 


rv = 4, we find E = 0.1011 x 10°. 
In general, E is 0 (hn?) as is the case with 


the discrete variable method. 


Error Analysis of the Method of Characteristics 


As discussed in Section C.1 of Appendix C, 
the trapezoid rule is used to integrate the charact- 
eristic equations. Consequently, the error, e, in 
the numerical computation is that incurred in the 


use of the trapezoid rule. 


Thus e is one where h = Ar = At. 
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Algorithm for Integrating Governing Equations for 
the Thin Compressible Shell 


The second order non-linear ordinary differ- 
ential equations governing the motion of the compres- 
sible thin shell are equations (3.30) and (3.31) 


written as follows. 


q = -2 [S - n26)26-2 | (A.13) 
h 
2 2 

UAB Mr 5 

h=qt-20- F) (A.14) 


To develop an algorithm for the numerical 


solution of these equations, let t, = kat, k = 0,1,2.. 


k 


represent the time after k equal intervals of At. 


v | 
Also, fom At) > 70, let hy = A(t) he <1 Bh (tidy 
hye = X (t,) and Oy = S(t), oe MER Ce Son ivbeael akong 
let hy be a function 
ye = A OAyey! 6p 4) (At LS) 
of A(ty_4) and 6 (t)_4)- 


Then under a known internal step input pres- 
sure g, there corresponds a by for each Ay,at any 
time t,: This by is found by solving the transcen- 
dental equation (A.13). The values of by and he 


are then used to integrate equation (A.14) by some 


Suitable stepwise procedure, for example the Pde 
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order Runge-Kutta process or the method of discrete 


variables. 


Equation (A.13) is best solved by the Newton- 


Raphson method. This, along with the discrete mech- 


anics procedure described in Section (A.1) is the 


basis of the following algorithm. 


Cisin) 


(iii) 


(iv) 


(v) 


(vi) 


(va) 


Cie) 


Read ei eal aly Apey = 0, and Ouay 


Spy is the initial guessed value for the 


= 1, where 


Newton-Raphson operation. 


Find the correct 6 Eromoa(Ae 3) sE0twathin 


Ket 


a specified error tolerance by the Newton- 


Raphson method. 


Use hyn] and Spy ey she bats! hyead from equation 
(A.15) 
ser Apel hye : 


Find an approximate A, using equation (A.2) 


k 


and from this, find an approximate i, using 


k 
equation (A.1). 


Substitute the approximate hye in equation 


(A.13) to find an approximate bye then, from 
(A.3), find approximate hie and denote this 
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(ix) 


(x) 


Ceaka) 


(xii) 


Pio 6 


Repeat steps (vi) to (vii) unt Ay is 
found within a specified error tolerance. 


Substitute the final he into equation (A. 2) 


t 
cOMEindmcthe.tinal hit ands this ine (Ane )eco 
find the final hee 

i] i] 
Transfer he ayes) Apeq! hy into heel ’ 


6, into by y LOe proceed to) the next. time 


and 


interval tears 


Repeat steps (iii) to (xi) for the interval, . 


le 


The flow chart for this algorithm is shown in 
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APPENDIX B 


CONDITIONS FOR A PERIODIC MOTION 


Periodicity of Motion 


Bel sieethe necessary condition 


A necessary condition for a periodic motion 
under a step function application of pressure, Q, 
is thatvan equilibriumsstace, o\ (ce), must) exisitor 
this pressure. The conditions for the existence of 
this equilibrium stretch is examined here for a 
thin shell for which \(e) refers to the mean equili- 
brium stretch. In what follows, we refer to the 
Mooney-Rivlin material as material I, the material 
with the logarithmic strain energy function of 
eguation (2.7) as material II, and the Blatz-—Ko 
material as material III. 

To investigate the existence of an equilibrium 
stretch for a given step input pressure, we examine 
the static curves, sketched in Figure 1, for the 
materials) l, Ll,) anG LiL. 

For material I, there are three distinct 
types of static curves marked as (i), (11) and) (1172) 
in Figure la. Curve (i) is the static curve for the 


neo-Hookean solid, (for which c = 1), curve (11) has 
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one maximum followed by a minimum point; this shape 


is obtained for all values of c for which 


re) xel Se Al (By 1) 


Cisne 
where Cre is the value of c at which the maximum 
and minimum points merge into an inflexion point de- 
noting a transition from curve (ii) to curve (iii). 
Curve (iii) shows Q increasing monotonically with ); 


this curve is obtained for values of c for which 


< < : 
0 c Cee (B.2) 


For the cases for which, 0 <cé<e there 


eran 
exists an equilibrium stretch for any value of the 
step input pressure Q, consequently, the motion re- 
sulting from any step input application of pressure 


is periodic. When, however, c <i Cos, there 


ebente 
could be more than one equilibrium stretch for a 
given applied pressure Q. Consequently, as discussed 
in Section (B.3), there exists a maximum value, Onur 
of the applied pressure Q, beyond which a periodic | 
motion is not possible. 

For materials II and III, the static curves 
are similar to the one for the neo-Hookean material 


and so, as in the latter case, there is a maximum 


value, Ores! of the applied step input pressure Q, 
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beyond which a periodic motion is not possible. The 


value of One is calculable as discussed in Section 


(Braye 


B.1.2 The sufficient condition 


The sufficient condition for a periodic motion 
i] 
is that the expression obtained for i - as in equa- 
tions (3.23) and (3.35), which can be represented 
24 


as a function, A = p(rA), of A, must be positive. 


Then and only then will the resulting motion be per- 
iodic, with two velocities, which are continuous 


functions of A thus 


A =f Yp(A) . (B. 3) 


This condition imposes some restrictions on 
the form of the strain energy function W, as shown 
by ZHONG-HENG and SOLECKI [11] for an incompressible 
material. 

Thus we have the following sufficient condi- 
tions to ensure a periodic motion of the incompres- 
Sible shell 


m 


W, OA) eM eas eos (B.4) 


and 


= 


W, (A) ~ NA as A} > 0 (B.5) 
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The Critical Value, c 


where M,N,m,n are positive constants such that m > 1 
andn?l, 
and where 


W, (A) = 5 wW(I,, I.) (B.6) 


rapes of c, for the Mooney- 


Rivlin Material 


The equation for the static inflation of the 
incompressible thin-walled shell of Mooney-Rivlin 


material, from equation (3.20), is 


2 
OuSIE 1 7 sae ah 
ce = 2c (i — os ap 21 Co) (r ey 
(B.7) 
Ors Lee Sia ee Sie 


The turning points of the Q/A curve, shown in Fig. 7, 


occur lat values Of 4 for which 


GQ = 9 


da if 
that is, at the roots of the equation 


92 - 2e(- H+ Fy + 2 - 1 + 3) = 0 
d d 


r 


Qu 
>hO 


(B. 8) 
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Putting 


and 


equation (B.8) becomes: 


i ee LS eo (B.9) 


By Descarte's rule of signs, this equation cannot 
have more than two real positive roots, so that there 
are at most two positive values of )\ for which Q has 
a stationary value. 

Furthermore, from «§(B.8), it follows that if 


c and (l-c) are both positive, that is, if we have 


(0) 4) fonsccowal (B.10) 
then 


se Soil Meeps uSeak 2G) ere SSG 


this implies that there exists some values of A’ > 1 
for which 


GQ < 


dr 


that is, for which from equation (B.9), we have 


cay Ay, 
f(x) = <1 Soar PTY. (Bele) 
x + 5x 


cepmoved (8.8) soltanpe 


: 2 
(@. a) v = - aye + 2 © BT 
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This is possible if and only if y is less than the 


maximum value which f(x) can attain. 


Solving the equation 


ohare 
dx 


= 0 
we find that the maximum value of f(x) occurs when 
x = 3.3883, that is, when A = 1.8407. Consequently, 


the critical value, c OL C11 sO. 8234. 


elionl ie 


The Maximum Step Input Pressure, Q , for a Periodic 


max 


Motion. 


Since for a periodic motion, the phase trajec- 


tory closes, the maximum value, Oe of the applied 


pressure Q, for a periodic motion is the value of Q 
for which the trajectory just fails to close as 
depicted in Fig. 2(5). | This limiting case occurs 
when an equilibrium stretch, Mle) o, coincides with 


a root, A of the energy equation 


a’ 
tue Verb epee 


Thus for the Mooney-Rivlin material, by setting 


A = 0 
es} ake) 
and ie 2 


in equations (3.20) and (3.23) respectively, we have 


at OF Bose al the pair of simultaneous equations; 
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2 
eto ae. eB Liita A Bete 
0 = 24 - 26( <3) AG ien =) 
(B.13) 
pes har ear rn egg shee eaves oe ee 
oe oT! 4 2,2 2 
(B.14) 


By eliminating Q between (B.13) and (B.14) 
we get a polynomial in A, the lowest positive root, 


Nae of which, when substituted in (B.13) gives 


Quan = 2e(eG - SS) + (LD - c) AQ - 
d 


max 
d 


> 
ale 


(BaL5) 
The values computed for Om for different 
cases of each of materials I, II and III are shown 
on the respective phase plane diagrams. 
Thea peciodsoteoscillation, T, is) finite ar 0 


is less than Os arg and is given by 
r 
T = ‘ an (B.16) 
CN 
The numerical method for evaluating this in- 


tegral is described in Program 5a. 
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APPENDIX C 


NOTES ON THE METHOD OF CHARACTERISTICS 


Integration of the Characteristic Equations 


Using the notations for the thick compressible 


Shell, let the functions 


a2 alee) 8 
eva 


ia = O,(r,T) and v = v(z,T) 


) 


(Cr) 


vary linearly in the small intervals between grid 
points. Then, referring to Fig. 8, the integral of 


the differential do for example, along Toss given 


[ do. (7, -(5,) (Or) 


and similar relations obtain for integrals along I 


by 


and II characteristics. 

Furthermore, by assuming the functions of 
equation (C.1) to be linear along the characteristics 
the trapezoid rule of integration becomes suitable 
for integration, from one grid point to the next, 
along these lines. For example, along the ine the 


following integration formula holds 
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(Ces3) 


For convenience of integration, we set 


Ay a) wt ee (C.4) 
Using equations (C.2) and (C.3) in’ (4.24) and 
(4.25), and dropping the bars on the dimensionless 
quantities, we get the integrated forms of the char- 


acteristic equations in non-dimensional form as 


follows. 
Along THe, 
) ( ) ) (2x 
a hae DOS Pees ihe 
ry ry 4 ak rj, r/, 
Oo Oo 
+ | | . (%2] Ar + B (:) + (z) INS Ge 
4 ik 4 ih 
(C75) 
along I , 
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i f D, pe 


along II, 
phe - Bc - (1 - eve, On ) 
= va we 
Siar ie + (i (2At) , 
(C2q7)) 
where 
sae 


and v is Poisson's ratio for infinitesimal strain 
from the undeformed state. 

To evaluate the stresses os and or and the 
Velocity Pveeat each gridipoint of thercharacteristic 
network of Fig. 4, we proceed as follows. In the 
eo. and’ Was aretaligizeror 


Is) 


Consequently, the values of these quantities at 


undisturbed region ABC, a 


grid points on line AC are obtained directly from 


the discontinuity relations of equations (4.26) - 
(4.28). With these initial values, a step-by-step 
solution is started by solving equations (C.6) and 
(C.7) simultaneously for the half-mesh in the lower 
corner of the triangular region ACD. This gives 
Omeande v at the third) pointwroneLneseolanguluar 


) 
mesh and since Ce is already prescribed at this 
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point, we can now solve equations. (C.5).= (C.7) 
simultaneously for the full interior mesh next to 
this first half-mesh. This gives values of Ov 4 
and Vv at the fourth point of the grid, from which 
we can evaluate the next mesh in a similar way. 
This process is repeated until all the grid points 
in the triangular region ACD have been treated. 

CD is a path of wavefront propagation, there- 
fore, the strong discontinuities in Or F% and. Vv. 
given by equations (4.26) - (4.28) must be added to 
the values previously calculated for these quantities 
along this line in order to proceed to the next tri- 
angular region CDE. CD is now the initial line and 
the foregoing stepwise process is repeated, this 
time starting with the first outer half-mesh in the 
lower corner of CDE. 

Jumps are again added for points along DE 
before going into the next triangular region which 
then becomes evaluable. In this way we can cover 


as many triangular regions as we desire. 
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' 
C.2 Boundary Conditions and the Constants K and K 


The discontinuity relations across an rae 


wavefront are 


- & 
[ols “= = 
ie hs as 
Ue ae =o er 
ay Aie sk 
JW] i. Seer (C.9) 
and across an I wavefront, these are 
1 
ag 
Koes ab ey 
i] 
Ait Vee se 
[9] a -v x 
1 
LEK 
[v] = + be a (Ge 40) 


where the notations for the thick compressible shell 
have been used and K, K are constants to be determ- 
ined from the boundary conditions. 
These boundary conditions are 
eae isp cue fe ee Gene Be ily) 


and 
o =O at r=B (or r=-). (Gerla) 
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For the first I wave propagation line AC, of 
Fig. 5, we have, from the first of equations (C.9) 
andethe  Elrst boundary conditions ins (C411); athe 


equation 
0 =IG = -q = = (Cal) 


where on o-. are the values of o, ahead Ofeandsbe- 


aro 
hind the wavefront respectively. Equation (C.12), 


therefore, gives 
K = -qA (Gr13) 


and the discontinuity relations along AC become 


= fee ish 
eine a 
= Vv _ GA 
[o,] a Saw ( =) 
© Seri) ets 
[v] = we eee (Crna) 


At the endpoint C of line AC, the first of equations 


(Cela) sqives 


and in order to satisfy thé second boundary condition 
in equation (C.11), we should have, from the first 


of equation (C.10), 


147 


— we 7 ~ ~~ 7 
7 - 7 — a 


io weil pigs 
Tr) rie 


wet Ee ok A 


ad hia Yo Obit 0 to ‘eevee ERS 8 tea” erty 
62.0) aadndigt | eft ia oerean sir algae edd Bake 
esvin ,exole weds 


ci.9) ad 


' ee "ad 
“. 


sacnnel tal Goole wiokieten —eeeemenekie aa tem 
Pe) = tot 

CP Ay = Tye 

(L.o) a ee | 


envl¥aves Ao dauct add , WA aatl do) % selodiine ene 2A 


eoidibage veh sn bags, aidoyTs Ldeevoy aaiere eb bas 
vael agit: harry vedealudds 2s 00) o) cel peopewat 
| (ah) neatiepe 20 


eh 


148 


This inplies thatgetonethe first I wave propagation 


Lane scp, 
ea (C15) 


and the discontinuity relations along this line 


become 
[o,}) = & 
yl ee eS 
Iv} = +5 (8). (c.16) 


Again, the first of equations (C.16) gives, 


for the endpoint D of CD, the discontinuity in o, as 


Consequently, to satisfy the first boundary condition 
in equation (C.11), we must have, for the value of 
K along the second ii propagation line DE, the ex- 


pression 
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K = -qA (8) 
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It follows, therefore, that for succeeding I propa- 


gation lines, the values of K are given by 
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K = -qA, +qA, -qA, +qA, ... (ents) 


with alternating signs, and for the I. propagation 


lines. 


Ket yrqht=qh;, +GA;7v-dqA, ... (C.20) 


also with signs alternating. 
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APPENDIX D 
PROOFS OF THEOREMS 


Proof of the Theorem Relating Surface Tractions 


To the Average Hydrostatic Stress 


The theorem leads to the equation 


| Yi BOPP ren tele, i | or * r aV = 3Vo (D.1) 
S Vid Wren VG 


~ ~ 


where T is the surface traction, r the position 
vector, V the volume, S the surface and o the average 
hydrostatic stress. 

Denote the left hand side of equation (D.1) 


by? =L;,-so that 
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and put) f= pr = Then, in itenson notation, equation 
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Applying the divergence theorem, this equation be- 
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By Cauchy's first law of motion, when there 


(Da) 


are no body forces present, we have 


It follows, therefore, from this and equation (D.3) 


that 


or 


as stated above. 


The Characteristic Curves 


Let the characteristic curve in the (r-t) 


plane be defined by the parametric equations 
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Then the dependent variables J OG and v have the 


following total differentials 
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thes partial derivatives, on Or O% eT Clamny, 
become indeterminate when the determinant of the 
coefficients of the left hand sides of the six 


equations (4.18) - (4.20) and (D.5) is zero. This 


condition is expressed by the equation 
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ac = +c and dr = 0 


where c is the velocity of propagation of dilata- 


tional spherical waves. 
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